Atom-light crystallization of BECs in multimode cavities: Nonequilibrium classical 
and quantum phase transitions, emergent lattices, supersolidity, and frustration 
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The self-organization of a Bose-Einstein condensate in a transversely pumped optical cavity is a 
process akin to crystallization: when pumped by a laser of sufficient intensity, the coupled matter 
and light fields evolve, spontaneously, into a spatially modulated pattern, or crystal, whose lattice 
structure is dictated by the geometry of the cavity. In cavities having multiple degenerate modes, the 
quasi-continuum of possible lattice arrangements, and the continuous symmetry breaking associated 
with the adoption of a particular lattice arrangement, give rise to phenomena such as phonons, 
defects, and frustration, which have hitherto been unexplored in ultracold atomic settings involving 
neutral atoms. The present work develops a nonequilibrium field-theoretic approach to explore 
the self-organization of a BEC in a pumped, lossy optical cavity. We find that the transition is 
well described, in the regime of primary interest, by an effective equilibrium theory. At nonzero 
temperatures, the self-organization occurs via a fiuctuation-driven first-order phase transition of 
the Brazovskii class; this transition persists to zero temperature, and crosses over into a quantum 
phase transition of a new universality class. We make further use of our field-theoretic description 
to investigate the role of nonequilibrium fiuctuations on the self-organization transition, as well 
as to explore the nucleation of ordered-phase droplets, the nature and energetics of topological 
defects, supersolidity in the ordered phase, and the possibility of frustration controlled by the cavity 
geometry. In addition, we discuss the range of experimental parameters for which we expect the 
phenomena described here to be observable, along with possible schemes for detecting ordering and 
fluctuations via either atomic correlations or the correlations of the light emitted from the cavity. 



I. INTRODUCTION 



Over the course of the past fifteen years, many phe- 
nomena of long-standing interest in condensed-matter 
physics have been realized in ultracold atomic set- 
tings IT]. Such realizations are considerably different 
from condensed-matter systems: in particular, ultra- 
cold atomic systems are highly controllable — i.e., they 
are isolated from the environment and are governed by 
thoroughly understood microscopic Hamiltonians — and 
tunable — i.e., the interaction strength, lattice depth, etc. 
are governed by quantities such as laser intensities, which 
are easy to alter. The high degree of control and tun- 
ability has made it possible both to explore emergent 
phenomena in a simpler setting than is typical in con- 
densed matter and to address hitherto experimentally 
inaccessible questions, such as the dynamics of ordering 
in systems that are quenched past a quantum critical 
point [5]. To date, most ultracold atomic realizations 
have focused on simulating the physics of electrons prop- 
agating through static lattices (via, e.g., realizations of 
the Hubbard model [3, ) or on constructing novel quantum 
fluids (e.g., To nks- Girardeau gases j^ or unitary Fermi 
gases [S]). Areas of condensed matter such as soft matter, 
supersolidity "5] , and glassiness — which involve emergent, 
compliant lattices capable of exhibiting dynamics, de- 
fects, melting etc. — have proven inaccessible to ultracold 
atomic physics because the lasers that create the lattice 
potentials in typical experiments are essentially insensi- 
tive to the atomic motion in those potentials. Aspects 
of such condensed matter phenomena remain unsettled 



in their traditional settings (e.g., the dynamics of glassy 
media and of supersolids), and therefore ultracold atomic 
realizations of them are especially desirable. 

A possible approach to realizing phenomena dependent 
on the emergent, compliant character of the lattice is 
to have the atoms interact with a potential created by 
dynamical, responsive quantum light, instead of static 
lasers. Exploring precisely such interactions has been 
the central theme of cavity QED [7] . Traditionally, cav- 
ity QED has aimed to realize systems involving a single 
atom coupled to a single mode of the electromagnetic 
field; however, the physics of many atoms coupled to one 
(or more) electromagnetic modes, i.e., many-body cavity 
QED, has also been studied extensively in recent work [S]- 
I16j . In particular, it was predicted in in Refs. ^."O] that a 
cloud of atoms confined in an optical cavity would exhibit 
collective effects such as self-organization; these effects 
were subsequently observed 110] . The atoms considered 
in these works were essentially classical, thermal parti- 
cles: however, the dynamics of many quantum atoms 
(e.g., a Bose-Einstein condensate or a Mott insulator) 
confined in a cavity has since been explored both theo- 
retically [I2l - [T4] and experimentally [HI [16]. Much of the 
work in this area, to date, has explored the novel implica- 
tions of the atom-cavity coupling for standard ultracold- 
atomic phenomena such as the superfluid-Mott insula- 
tor transition |T3l [14] or the collective excitations of a 
Bose-Einstein condensate [T21[IS]. The objective of the 
present work is to suggest that a quite different class of 
condensed-matter problems, involving the emergence and 
dynamics of spatially ordered states, can be realized and 
explored using ultracold atoms confined in optical cav- 



ities. Elements of this work were reported in Ref. [T7] : 
related issues, involving the simulation of phonons in op- 
tical lattices, were previously raised in Ref. [T5] . 

A central idea in the extant literature as well as 
the present work is the idea of cavity-induced self- 
organization [HI [9], which may be explained as follows: 
Consider N two-level atoms in a single-mode optical cav- 
ity, interacting with the cavity mode and a pump laser 
oriented transverse to the cavity axis (see Fig. IT]). The 
atoms coherently scatter light between the pump and 
cavity modes. Atoms arranged at every other antinode 
of the cavity field (i.e., one cavity-mode wavelength A 
apart) emit in phase; therefore, A-period fluctuations of 
the atomic density increase the number of photons in the 
cavity, thus drawing the atoms into A-spaced wells at ei- 
ther the even or odd antinodes. This leads to greater 
constructive interference in the emitted light, stronger 
atomic trapping, and so on. The system reaches a spa- 
tially modulated steady state when the energetic gain 
from the atom-light interaction is balanced by the cost, 
in kinetic energy or repulsive interactions, of confining 
the atoms to either the even or the odd sites of the emer- 
gent lattice. 

Although self-organization in a single-mode cavity re- 
sults in the spontaneous breaking of a discrete symme- 
try, the locations of the antinodes are not themselves 
emergent, but are fixed by the cavity geometry. In other 
words, the non-crystalline state does not possess contin- 
uous translational invariance; thus self-organization re- 
sembles, e.g., a phase transition between crystal struc- 
tures, rather than true crystallization. With multi- 
mode cavities, by contrast, self-organization results in the 
breaking of continuous symmetries, both in the collective 
choice of which mode(s) to populate with photons and in 
the choice of relative phases between the modes. For ex- 
ample, in the case of the ring cavity (which consists of 
two counter-propagating traveling- wave modes [H]), the 
atoms must collectively choose the (continuous) relative 
phase between the two counter-propagating modes, thus 
setting the location of the antinodes of the cavity field. 
As with real solid-state crystallization, the breaking of 
this continuous symmetry induces rigidity with respect to 
lattice deformations. With larger families of modes, one 
can envision realizing such characteristically crystalline 
notions as dislocations and geometrical frustration. 

In a previous article [T7] , we developed a field-theoretic 
description of the interacting many-atom, many-mode 
system, which we applied to the case of a transversely 
pumped concentric cavity. We observed: (i) that a quasi- 
two-dimensional atomic cloud undergoes a weakly first- 
order transition into a spatially ordered state, and that 
this transition becomes a quantum phase transition at 
T = 0; and (ii) that for the case of a strongly layered 
three-dimensional cloud, inter-layer frustration precludes 
global ordering, and the system instead breaks up into in- 
homogeneous, static domains. In addition, we suggested 
that both fluctuation phenomena and signatures of su- 
persolidity may be observed via the spatial and temporal 



correlations of the light emitted from the cavity. The 
formalism developed in Ref. }17j assumed, however, that 
because the flux of energy through the atom-cavity sys- 
tem is negligible in the regime of interest, one could adopt 
a quasi-equilibrium description of this phase transition. 
One of our objectives in the present work is to justify this 
assumption within a general nonequilibrium formalism; 
furthermore, we compute the leading corrections to the 
effective theory of Ref. [TTj that arise because of nonequi- 
librium effects. 

In the present work, our approach towards meeting 
these objectives is as follows: (i) we develop a fully 
nonequilibrium, field-theoretic description of the atom- 
cavity system, using the Schwinger-Keldysh functional- 
integral formalism [19l[20]; (ii) we show that the nonequi- 
librium description can be reduced, in the regime in 
which the cavity's photon decay time is longer than the 
timescales for atomic motion, to an effective equilibrium 
description; (iii) we analyze this effective equilibrium the- 
ory using diagrammatic and renormalization-group tech- 
niques to establish the nature of the self-organization 
transition, in the specific case of a concentric cavity; 
and (iv) we reintroduce the nonequilibrium effects, due 
to the leakage of photons out of the cavity, using per- 
turbation theory, and account for their effects on criti- 
cal behavior near the self-organization transition. Our 
analysis of the equilibrium theory extends that used in 
our previous work |17| via an adaptation of Shankar's 
renormalization-group treatment of the Fermi liquid |21j , 
as well as an analogy with the 0(p)-invariant vector 
model of magnetism, to establish that for an interacting 
Bose-Einstein condensate the self-organization transition 
is always discontinuous in the concentric cavity. Further- 
more, we show that the chief consequence of nonequi- 
librium effects is to decohere quantum correlations on 
a timescale related to the linewidth of the cavity; thus, 
the self-organization transition is always classical on the 
longest timescales. (For sufficiently high-finesse cavities, 
there should, however, be an extended crossover regime of 
timescales for which the phase transition appears quan- 
tum.) 

After presenting our analysis of the phase transition 
itself, we turn to the properties of the ordered (i.e., 
self-organized) state. We show that the ordered state 
has low-lying excitations, associated with the continu- 
ous symmetry-breaking, that resemble the excitations of 
smectic liquid crystals. We also expand on our previ- 
ous observation [17] that, for a Bose-condensed atomic 
cloud, the ordered state would be a "supersolid" (or, 
more accurately, a "super-smectic" ) in that it simulta- 
neously possesses emergent (liquid-)crystalline and su- 
perfluid order. The properties and even the existence 
of supersolidity in ^He are much-discussed topics in the 
condensed-matter literature [51 JHHIS]; the appeal of 
ultracold-atomic realizations is that one can explore the 
characteristic phenomenology of supersolids in contexts 
where it is less challenging to establish that they do 
in fact possess supersolidity. In order to explore the 



relevant phenomenology, it is necessary that the solid- 
ity be associated with a broken continuous spatial sym- 
metry; for this purpose a continuum supersolid such as 
the one explored in the present work should serve as a 
more suitable setting than would, e.g., the lattice super- 
solids proposed in Refs. [SSI EZ]- The present scheme 
has the added advantage of being more readily realiz- 
able. The self-organization of a BEC in a cavity was, in 
fact, recently demonstrated experimentally by Baumann 
et al. |16j for the case of a single-mode cavity; generaliz- 
ing this experiment to the multiniode case, in which one 
has continuous symmetry-breaking, should be technically 
straightforward. Continuing with the theme of superso- 
lidity, we propose a scheme for detecting such as state, 
and develop a schematic phase diagram for the system 
(Fig. 11), which exhibits three phases: the supersolid. 



the normal solid, and the uniform superfluid. The phase 
diagram for a multimode cavity differs from that for a 
single-mode cavity in that the multimode case features a 
direct transition from the uniform superfluid to the nor- 
mal solid, whereas in the single-mode case there is always 
a supersolid regime separating the uniform superfluid and 
the normal solid. 

This paper is organized as follows. In Sec. [ll] we de- 
scribe the microscopic model of the atom-cavity system 
that is used in the rest of the paper, and in Sec.|in]we dis- 
cuss the qualitative behavior one might expect from this 
model. In the next three sections we construct and an- 
alyze a nonequilibrium field-theoretic formulation of this 
model: Sec. llVl introduces the relevant field-theoretic for- 
malism. Sec. [V| applies this formalism to derive an atoms- 



only action, and Sec. VI describes the quasi-equilibrium 
limit of the atoms-only action. In Sec. |VII| we derive 
an effective Ginzburg-Landau free energy, valid near the 
phase transition, which realizes a version of Brazovskii's 
model ,28j of ordering at a finite wavelength, and an- 
alyze the effects of fluctuations on the self-organization 
transition in both the classical (T > 0) and quantum 
(T = 0) cases. In Sec. VIII| we turn to the effects of 
departures from equilibrium, both on the fluctuations 
near the transition and (following the work of Hohen- 
berg and Swift [H]) on the nucleation of ordered states. 
In the next two sections we focus the properties of the 
ordered state: Sec. |IX| reviews the properties and ele- 
mentary excitations of the ordered state, and Sec. |X] dis- 
cusses the supersoHd aspects of the ordered state. In 



Sec. XI we discuss the experimental feasibility of the phe- 
nomena that we are investigating, showing that most of 
them should be readily detectable in the laboratory. In 
Sec. |XII| we briefly consider the case of strongly layered 
three-dimensional systems, which have the feature that 
geometrical factors tend to frustrate the development of 
globally coherent long-range order. Finally, in Sec. |XIII| 
we summarize the results of this work, and discuss its 
relationship with other problems involving phase transi- 
tions and related collective effects in atom-light systems. 
Various supplementary issues are addressed in four ap- 
pendices. 



II. MODEL 

The system analyzed in this paper consists of N two- 
level atoms confined in a multimode optical cavity, to- 
gether with the electromagnetic modes of the cavity. The 
system is pumped by an external laser, which is oriented 
transverse to the cavity axis, as shown in Fig. IT] In addi- 
tion, the system is coupled to a set of extracavity electro- 
magnetic modes, which constitute a bath for the system. 
The complete Hamiltonian H governing the system and 
bath comprises three elements — viz., -ffat, the atoms-only 
Hamiltonian; i?cmi the light-only Hamiltonian; and -ffint, 
the atom-light interaction Hamiltonian — which we dis- 
cuss in detail in the rest of this section. 



1. Two-level atoms 



The atoms are described by the Hamiltonian 
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where N is the number of atoms (indexed by n = 
1, . . . , iV), each of which has mass M; the position and 
momentum of atom n are, respectively, x„ and p„. The 
operator a^ denotes the Pauli operator, which acts on the 
internal state of the atoms, which, for the sake of sim- 
plicity, we take to be two-state atoms, with fiwA being 
the energy splitting between the ground (g) and excited 
(e) states. (Our conclusions do not hinge in any essen- 
tial way on this restriction to two states.) We model the 
interaction between atoms via a repulsive contact po- 
tential, which is parametrized in terms of the (positive) 
parameter U . To ensure the correct handling of the Bose- 
Einstein quantum statistics of the atoms, we employ the 
framework of second quantization. Thus, H^^ becomes 

H,t=-^Jd''x{^l{^)VH^i^) + vl.t(x)v2*^(x)} (2) 

+hu;A /"d'*a;*^(x)*,(x) 
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in which the pair of bosonic field operators, ^l'g(x) and 
\I'e(x), respectively represent the ground- and excited- 
state atoms. In the last line of Eq. (l2]), all the coordinate 
indices are x, as the interaction is taken to be local in 
space. 

Next, we assume that the density of excited-state 
atoms is sufficiently low that collisions between such 
atoms can be neglected. This assumption necessarily 
holds in the low-temperature regime, which is the regime 
of primary interest to us, because the rate of spontaneous 
decay (which is proportional to the number of excited 
atoms) must be kept low in order to avoid the heating 



associated with it. With this regime in mind, we approx- 
imate the interaction term as 



U 
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(3) 



In what follows we shall replace U by the frequency U = 
U /h, in order that its magnitude be expressed in the 
same units as those conventionally used to describe the 
atom-cavity coupling. 



2. Electromagnetic modes 

Optical cavities typically consist of two or more mir- 
rors that support localized modes of the electromagnetic 
field between them [3D]. The modes might be stand- 
ing waves, as in the single-mode cavity, or the traveling 
waves that one can have in the ring cavity. In general, 
we shall be concerned with transverse electromagnetic 
(TEM) modes. The vector character of the electromag- 
netic field can be absorbed into the effective atom-mode 
coupling 0; hence, such modes can effectively be de- 
scribed in terms of harmonic solutions to the scalar wave 
equation [30] . Furthermore, as the mirrors are not per- 
fectly reflective, these modes are in fact weakly coupled 
to the extracavity modes, and thus cavity-mode pho- 
tons tend to "leak" out of the cavity at some nonzero 
rate k (which sets the intrinsic linewidth of the cavity) . 
These considerations lead us to model the electromag- 
netic field — both intracavity and extracavity — in terms 
of the Hamiltonian 
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where Oq and A^ respectively represent the intracavity 
and extracavity photons, uja and w^ are the intracav- 
ity and extracavity mode frequencies, and n^^e describes 
the coupling between the intracavity mode a and the 
extracavity mode e. We assume that the intracavity- 
extracavity coupling is weak enough (i.e., the cavity is of 
sufficiently high finesse) that it is meaningful to separate 
the modes into intracavity and extracavity ones. 

The modes of a generic standing-wave cavity are not 
frequency-degenerate: the typical frequency spacing, or 
longitudinal "free spectral range," is of order c/ L, where 
L is the linear dimension of the cavity: e.g., for a 1 cm 
cavity, the free spectral range is about 15 GHz. For 
such a cavity, the frequency spacing between higher- 
order transverse modes is an appreciable fraction of this 
number; thus there are no degenerate modes. However, 
certain specific cavity geometries do support degenerate 
modes. The simplest of these is the ring cavity, which 
is a three-mirror arrangement that supports two coun- 
terpropagating traveling- wave modes (see Fig. flk) . Even 



larger degeneracies are possible, e.g., in confocal or con- 
centric cavities. Let us label the cavity-mode structure 
by the integers {I, m, n), where n is the number of nodes 
along the cavity's axial direction, and (/, m) are the num- 
bers of nodes along the transverse directions; see Fig. [T] 
(The corresponding mode functions are approximately 
sinusoidal in the axial direction, and Hermite-Gaussian 
(or Laguerre-Gaussian) in the transverse directions, al- 
though this approximation breaks down in the limit of a 
concentric cavity.) In the confocal cavity, the condition 
for frequency degeneracy is that n + {{I + m)/2) = M 
(see, e.g., Ref. [3D]) for some fixed integer M; in the con- 
centric cavity, the condition becomes I + m + n = M . In 
principle, these conditions imply that there are of order 
A/^ degenerate modes, where M is roughly the number of 
optical wavelengths across the cavity, commonly 10** or 
more. In practice, higher-order modes are increasingly 
lossy, because their profiles (i.e., spot sizes) are larger, 
and therefore more of their light leaks out of the sides of 
the cavity mirrors, which typically occupy only a modest 
amount of solid angle. We approximately account for this 
effect by assuming that Mq of the modes have a common 
loss rate k, and that the other modes are perfectly lossy 
(i.e., for them k = oo). 



3. Atom-ligfit interactions 



In the dipole and rotating- wave approximations [7] , the 
atom-light coupling has the generic form zct_ ajj5„(x) — 
h.c, where the a operators raise or lower the atomic in- 
ternal state; a„ is either an intracavity or extracavity 
mode; and 5n(x) = (7S„(x), in which g is the coupling 
between the atom and the mode, and S„ is an appro- 
priately normalized mode function. Note that g is pro- 
portional to 1/VV, where V is the system volume |7j; 
hence g'^N stays finite in the thermodynamic limit. In 
the special case of the pump laser mode, we shall treat 
the corresponding a„ as a classical variable, so that the 
atom-laser coupling takes the form ri(x, t)a-+]i.c., where 
fi(x) = f7(x)e-'"^*, with fl(x) being the local Rabi fre- 
quency (which is proportional to the local amplitude) of 
the laser. In our analysis, we shall largely neglect the 
term that governs spontaneous atomic emission into ex- 
tracavity modes (i.e., spontaneous decay); this term is 
proportional to the linewidth of the transition, which is 
denoted as 7. The reason we can neglect spontaneous de- 
cay processes is not that such processes are always weak. 
Rather, it is that their effects amount to the heating up 
of the sample via the random impulses they give to the 
atoms, and therefore the timescale for spontaneous decay 
[which is given by 1/Rj = A\/{'-fD,'^)] acts as an upper 
limit on the duration of an experiment: quantum dynam- 
ics on timescales slower than R^ is likely to be washed 
out as a result of spontaneous decay. Neglecting, then, 
spontaneous decay processes, the Hamiltonian governing 
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FIG. 1. (a) The transversely pumped, quasi-two-dimensional geometry primarily discussed in this paper, (b) Ring cavity 
geometry. The pump laser beam is perpendicular to the plane defined by the three mirrors, as indicated in the figure, (c) 
Schematic representation of a concentric cavity, showing the partial rotational symmetry that such a cavity inherits from the 
sphere of which both cavity mirrors are arcs, (d) Three-dimensional view of a representative mode function for the concentric 
cavity. This mode function is labeled by {l,m,n) = (2,1,5), or alternatively by TEM21. The three numbers enumerate the 
nodes (one fewer than the number of lobes) in the pump (z), angular, and radial directions, respectively. The axial mode index 
n is fixed by the requirement that I -\- m -\- nhe constant for a family of degenerate modes, and can therefore be suppressed, 
(e) The intensity profile of the representative mode TEM21 at one of the cavity's end mirrors, (f) The intensity profile of the 
mode TEM21 in the equatorial (i.e., z = Q) plane of the cavity. 



the atom-light interactions is given by 
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III. PHYSICAL EXPECTATIONS 

A. Single-mode case: semiclassical picture 

As discussed in the Introduction, the atom-cavity sys- 
tem is unstable towards crystaUization when transversely 
pumped at sufficient intensity. A more quantitative, al- 
beit semiclassical, picture of the crystallization is as fol- 
lows. Assume that the atoms are pumped by a pair of 
counterpropagating lasers perpendicular to the equato- 
rial plane of the cavity (see Fig. flk) , so that the electric 
field due to the lasers is given by El cos kzy, and that 
the field in the cavity mode, which is a standing wave, is 
given by Ec cos kxy. The polarizations of the two modes, 
being parallel, can be neglected for the present purposes. 



Atoms in the plane z = Q are subject to an effective elec- 
tric field of intensity E\ -\- 2ElEc cos kx + E"^ cos^ kx. 
The first term is spatially constant; of the spatially vary- 
ing terms, for small Ec/E^ the second is the dominant 
one. The potential energy of the (high-field-seeking) 
atoms in this field can therefore be taken to be 



£ ^ —Ec'El / dxn{x) cos kx. 



(6) 



Furthermore, the magnitude of Ec is set by the atomic 
distribution, i.e., 
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where 
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(j) is the pump laser-cavity phase differ- 
From Eqs. (|6| and ^ it follows that £ ~ 
—rij,. When £ exceeds the 
kinetic-energy cost of self-organizing, the atoms undergo 
self-organization. As £ depends only on nf., it is invari- 
ant under rtfc — )• — n^, i.e., under translating the distri- 
bution of atoms through half a mode wavelength (e.g., 
from even to odd antinodes) along with changing the 



sign of the cavity field. This is the even-odd symme- 
try (which is an inversion symmetry, provided the cavity 
has an even number of total wavelengths) that is broken 
when the atom-light system crystallizes. This crystalliza- 
tion transition can be thought of in one of two equivalent 
ways: either as a crystallization of the atoms or as a 
locking of the phases of the laser and the cavity modes. 
This phase-locking is associated with the condensation 
of photons into a cavity mode; analogously, conventional 
crystallization can be thought of as a condensation of 
phonons. 

Note that, for a single-mode cavity, the presence of 
the term in £ that goes as Eq cos^ kx does not quali- 
tatively affect the picture just outlined, unless Ec is of 
the same order of magnitude as Ej^. If Ec > Ej^, the 
energetics would still favor symmetry breaking; however, 
the dynamics of ordering would then be complicated by 
the fact that there are local electric-field minima at the 
"minority" antinodes — i.e., the ones at which the atomic 
density is supposed to be low — and these local minima 
can trap atoms, as discussed in Ref. [9]. We shall not 
discuss this regime further in the present work. 



B. Multimode case 

The question that arises when one attempts to extend 
the idea of atom-light self-organization to the setting of 
multimode cavities is this: Which mode(s) do the atoms 
crystallize into? In this section, we offer some heuristic 
general considerations aimed at addressing this issue. We 
shall return to this question in the specific context of the 
concentric cavity, once we have developed the relevant 
field-theoretic techniques. 



1. Traveling waves 



The argument given in Sec. |III A| on self-organization 
for a single-mode cavity proceeds similarly for the case 
of multimode cavities, except that the individual mode- 
function coskx must be replaced by the (as yet unspec- 
ified) set of cavity mode functions {(^^(x)}. There are, 
however, two crucial differences. The first is pertinent 
whenever the cavity supports traveling- wave modes (e.g., 
the ring cavity). In this case, the potential energy is given 

by 



d'':c5,,(x)n(x) / ^^^^(x') n(x'), (8) 



and the dynamics of each mode is coupled to that of 
its partner under time-reversal. This has an important 
consequence, which is easiest to illustrate in the case of 
the ring cavity. Here, ga = expikx, and £ is invariant 
under the transformation n{x) — > n{x + e), which in- 
volves shifting the atomic distribution along the cavity 
axis (and adjusting the antinodes of the cavity mode ac- 
cordingly). Therefore, crystallization in multimode cavi- 



ties having traveling-wave modes necessarily involves the 
spontaneous breaking of a continuous translational sym- 
metry. 



2. Mode selection 

If two cavity mode functions ga (x) and ga' (x) are as- 
sociated with a pair of frequency-degenerate harmonic 
solutions of the wave equation for the same (homo- 
geneous) boundary conditions, any linear combination 
Cq (7q(x) -f Ca' ga'(x) is also a legitimate cavity mode. 
Therefore, it might seem that, in an iV-fold degenerate 
cavity, any normalized mode of the form X]a^a5a(x) 
would be an "equally good" arrangement for crystal- 
lization, i.e., there is an A'^-dimensional degenerate sub- 
space. This is not generally true, as a result of terms 
in the energy, omitted so far in the present section, that 
lift this degeneracy, such as the interatomic contact re- 
pulsion. Consider the extreme simplification involving 
two modes having the respective mode functions cos kx 
and cosky (as would arise, e.g., from two cavities, per- 
pendicular to one another and to the laser): any func- 
tion of the form {Ca,Ca') = CaCOs{kx) + Ca' cos{ky), 
with Ca + C^, = 1, is a legitimate mode function. 
If the atoms are self-organized in the state {Ca,Ca'), 
the expectation value of the atomic field is given by 
(ipi^X-)) ^ A + B [Ca cos{kx) + Ca' cos(fcy)], and the in- 
teraction energy, Eq. ([5|, goes as / d'^x |^(x)|^. This can 
readily be checked to be smallest when either Ca — or 
Ca' = 0, i.e., for a stripe-like arrangement along either 
the X axis or the y axis. 

A similar effect arises from the mode-mode scattering 



term (i.e., the E^ term). 



d'^xn{x.)ga{x.)ga'{^)- 



(9) 



In the two-mode example discussed in the previous para- 
graph, in which the modes are at right angles to one an- 
other, this term is essentially diagonal in the mode indices 
for either of the stripe-like states. Suppose, however, that 
the two cavities lie at a small angle 6 rather than at right 
angles to one another, so that the modes are cos(k • x) 
and cos(k' • x) with k « k'. In this case, atomic density 
fluctuations of wave-vector |k — k'| « |k|0 — which, for 
small 6, could be excited either thermally or quantum- 
mechanically — would suffice to mix the cavity modes. 
The effect of such mixing would be to lock the relative 
phases of the two modes. 



IV. FIELD-THEORETIC FORMULATION 

Our objective in this and the subsequent two sections. 
Sec. [V| and Sec.|VH is to construct a useful field-theoretic 
formulation that will enable us to explore the quantum 
statistical mechanics of correlated many-atom, many- 
photon systems in multimode cavities. Having done that. 



in Sec. VII we employ this formulation to address issues 



such as the emergence and nature of the spatial struc- 
ture and spatio-temporal atomic and photonic correla- 
tion properties of such systems, focusing on the vicinity 
of the transition to the self-organized state. The atom- 
cavity systems of interest here are neither closed nor in 
thermal equilibrium, because they are driven by an exter- 
nal pump laser (which adds energy) and leak photons into 
the continuum of modes that lie outside the cavity (hence 
losing energy); thus, even in its steady states there is a 
flux of energy through the system. For these reasons, the 
structure and correlation properties must be computed 
within a nonequilibrium formalism. The one we employ 
is the closed-time-path formalism, due to Schwinger [19] 
and Keldysh POJ , which enables the use of diagrammatic 
methods as well as renormalization-group techniques to 
analyze fluctuations. (For a discussion of the differences 
between the equilibrium and nonequilibrium formalisms, 
seeSec.|v!]andRef. [S]-) 

Although, as we have just discussed, a full analysis 
of the problem demands a nonequilibrium approach, we 
find that, for systems that are near the threshold for self- 
organization and in the dispersive regime (i.e., the pump 
laser is far-detuned from the atomic resonance), an ef- 
fective equilibrium description is valid, to a reasonable 
approximation. As we shall see, e.g. in Sec. |VIII[ the 
description of this regime can be improved, and other 
regimes (such as the strongly organized regime) can be 
analyzed, using the full machinery of the Schwinger- 
Keldysh nonequilibrium approach. 



their response and correlation functions. Formally, the 
task of computing these may be stated as follows: sup- 
pose that we know the state of the system in the infi- 
nite past, as described by its density matrix p{t — — oo), 
when it is taken to be isolated and noninteracting. The 
system is then coupled to an environment (or environ- 
ments) , which generically force the composite of system 
and environment (s) to be out of equilibrium, and the 
intra-system interactions are adiabatically switched on. 
The question then becomes: What are the expectation 
values and response and correlation functions of the var- 
ious system observables, once the system and environ- 
ment have relaxed to a steady state? 

Let us first consider the case of a single harmonic oscil- 
lator with Hamiltonian H — HutQ b^ b, i.e., a free bosonic 
degree of freedom having the characteristic frequency ujq. 
Suppose, as a simple example, that we are interested in 
the expectation value of some observable A(t) at time t, 
given that the system was at some time ti < t in a ther- 
mal state at temperature T, i.e., governed by the density 
matrix p{ti) oc exp(— /two^^fo/fc^T). Thus, we wish to 
compute the quantity 



(^(t))^Tr[p(t,)A(t)] 



(10) 



A. Schwinger-Keldysh functional integral 

The quantities of interest in quantum many-body dy- 
namics are the expectation values of observables and 



One can expand the trace in terms of bosonic coherent 
states [32] at the times ti and tf > t, thus arriving at the 
expression 



(11) 
Note that the primary motivation for inserting the complete set of states \wf) is to make the above expression more 
symmetric between initial and final times. The relevant matrix element of the initial density matrix is given by 
e:iq3{— /Shujow* Wi) , the overlap {wf\l\wf) is given by exp{—w%Wf), and the two other expressions, which are transition 
amplitudes, can be rewritten as coherent-state path integrals: 



{w,\wf){wf\A{t)\w,) = Z?(z+,z;)e^^[^+^<l / i?(z_, zl)A(z_, zl)e-''^[^-^^* 1, (12) 

J Z^{ti)—Wi J Z—{ti) — Wi 



in which the ± signs indicate whether the final state is the ket (-I-) or the bra (— ), and the action S is given by 



S[z±,zl.]= / dtz^{idt ~uJo)z±. 



(13) 



One can thus rewrite Eq. ( 10 ) as follows: 



(^W> = / d{wfW*fW^w*WiW*Wfw))e- 



\'W f\ —\'Wi\ —\wf\ —\'Wi\ —jShujQW^W^ —Wj.-W 



D{z+,zl) / Diz^,z*_)-—- 



If we omit the functional derivative S/6^-{t) from the 
right hand side of Eq. ( |14[ ), the remaining formula is 
the quantity commonly denoted as Z (by analogy with 
the partition function). Z is a generating functional for 
the correlation functions of the oscillator: it can be dif- 
ferentiated repeatedly with respect to either of the two 
source functions ^± (t) in order to generate all requisite 
correlation functions of A. By differentiating with re- 
spect to the -|- and — sources appropriately, one can 
compute expectation values that involve various order- 
ings of the operators, e.g., the time-ordered, retarded, 
and advanced Green functions. By contrast, in the zero- 
temperature and Matsubara nonzero temperature equi- 
librium formalisms, the only correlation functions that 
can be computed directly are the time-ordered ones; the 
(physically relevant) retarded Green functions are then 
to be inferred using identities that hold in equilibrium or 
at zero temperature. 



Because, in Eq. (14 1, the initial and final values of the 



paths are integrated over (with an exponential measure) , 
the path integral in Eq. ( |14[ ) is in effect an unconstrained 
path integral over the two sets of paths z±. Moreover, 
these paths are uncoupled from one another except at the 
two endpoints of the path integrals. For the oscillator in 
question one can write the action in the form 



S = 






S+- 






(15) 



in which integration is implied over time. We shall 
sometimes denote as S the matrix comprising the four 

block matrices S±±. The diagonal blocks 5*++ and S 

are given by Eq. (13); the off-diagonal blocks S^ and 



S (- are zero, except at the time-endpoints; they can- 
not, however, be neglected, because the presence of such 
off-diagonal terms changes the inverse of S, denoted G, 
which contains the two-point correlation functions and 
has the following structure: 

G+_(t,t')=nB('^o)e-'"°(*-*') 

= {z{t)z*{t')), (16a) 

G_+(<,t') = (ns(c^o) + l)e-^"°(*-*') 

= {z*{t)z{f)), (16b) 

G++{t, t') = <d{t - t')G+^ + e{t' - t)G^+, (16c) 

G__(i, t') = e{t - t')G-+ + e(i' - t)G+-, (16d) 

where ^^(wo) is the Bose-Einstein distribution function. 
Evidently, G_|__|_ and G are the time-ordered and anti- 
time-ordered Green functions. It is convenient for our 



e 



^-aifOi; 



iS[z+,zl]-iS[z^,z'_]+} dtS,-(t)A-(t)+^+(t)A+(t) 



,(14) 



r 



purposes to rotate z± into its "classical" and "quantum" 
components, defined as follows: Zc = [z^ + z_)/2, Zq = 
(z+ — z^)/2, in Eq. (15). This has the advantage of 



reducing the number of independent Green functions by 
one: 



where 



Gk = 
Gr = 
Ga = 



r^ _ I Gk Gii 

^-^ Ga 



«(t)z,(i'))=-*(2nB(^o) + l)e 
{z*it)z^it'))=-^eit-t')e-^-°('-''\ 
{z*{t)z^{t'))^ ze(t'-i)e— o(*-0. 



(17) 

''^oi*-^')^ (18a) 
(18b) 
(18c) 



Gji and Ga are the retarded and advanced Green 
functions — which describe the response of the system to 
an external perturbation — whereas Gk, the "Keldysh 
Green function',' depends on the system's initial den- 
sity matrix and its correlations. For a system in equi- 
librium with a bath, the correlations and response are 
related by the fluctuation-dissipation theorem; by con- 
trast, for an isolated, noninteracting (and hence nonequi- 
librating) system such as the harmonic oscillator, the 
two properties — response and correlation — are indepen- 
dent. In interacting systems that are away from equi- 
librium, there is in general a complicated interplay be- 
tween the correlations and the response; therefore, all the 
Green functions contain information about both correla- 
tions and response. The nature of this relation, however, 
varies from system to system, and hence the information 
contained in the three Green functions is not redundant. 



B. Application to atom-photon system 

The prescription for proceeding from the second- 
quantized Hamiltonian for a generic set of bosonic de- 
grees of freedom {0i(i)} to the appropriate coherent-state 
path integral (in real time) is as follows: 

-l-^M,|0,,+ (t)p-(-|-o-). (19) 

r 

In this expression, /i^ is the chemical potential for field 
1 — in the present case, the chemical potential for the pho- 
tons is zero, whereas that for the atoms determines the 



atomic density — and the symbol (+ ^ — ) indicates a 
formally identical set of terms in which + fields are re- 
placed by the corresponding — fields. This prescription 
can be carried out for each of the terms in the Hamilto- 
nian H introduced in Sec. [nj and generates all the ++ 

and components of the action. The off-diagonal (H — 

and — h) blocks in the bare (i.e., microscopic) theory de- 
pend on the appropriate Bose-Einstein distributions of 
the free photonic modes and the free atomic modes. Be- 
cause these blocks are infinitesimal (as they arise from the 
end-point couplings discussed in the previous section) , it 
is more useful to write down the relevant blocks in the 

inverse action, viz., the bare Green functions G^ and 

G 1-. For the cavity photon modes these have precisely 

the forms in Eqs. (16), i.e.. 



G+_(i,t') = «B(^c)e-*""^*-*'\ (20) 

G_+(i,t') = {riBiuc) + l)e— ^(*-*'), (21) 

whereas for bosonic atoms (in their internal ground state) 
and in a single-particle eigenstate of the kinetic energy 
having eigenvalue E, these have the form 



G+-{t,t') =nB(^-^)e-^(-^-'^'(*-*')/'^ (22) 

G_+(t,t') = {uBiE -^l) + i)e-«(s-M)(t-t')/^ (23) 

Other degrees of freedom can be treated similarly. 



CONSTRUCTING THE ATOMS-ONLY 
ACTION 



In this section we derive an effective action, involv- 
ing the ground-state atoms, which we shall use to deter- 
mine expectation values and correlators involving atoms 
and/or intracavity photons. This is accomplished by in- 
tegrating out all other degrees of freedom — a task that 
is straightforward, owing to the fact that they appear 
quadratically in the complete action. 



A. Eliminating the atomic excited state 



As we see from Eqs. p^, (pi), and (Is]), the complete action involves the atomic excited state in the following terms: 



dtd'^x^l^{:>i,t) (id. 



2M 



UJA 



*e,±(x,i)-Ki[^.9a(x)*g,±(x,t)*];^±(x,i)a„,±(i)-h.: 



-C/|vI/<,,±(x,i)nvI/,^±(x,t)|2 + 



(24) 



The functional integral over the quadratically occurring '^e,± hi Eq. (24) can be performed exactly. In the regime of 
interest, the atom-laser detuning HA a = huJA — fn^L is much greater than the energy scale associated with atomic 
motion; therefore, one can simplify matters by dropping the gradient term for ^e- (Put heuristically, excited-state 
atoms, being short-lived and massive, "decay" before they have time to move, so that the interactions they mediate 
are local in space and time.) Thus, one can integrate out the excited state, determining the necessary kernel via the 
standard technique of solving the classical equations of motion for ^e.± and 'i'], j. to obtain 



*e,±(x,t) 



.EagaWaa(t)^g,±(x,t) 



-iuj^t 



(25) 



to enable us later to exploit the 



where, for convenience, we have made the change of photon variables aa — >■ aa e 

approximate degeneracy of the laser and cavity modes. By inserting the classical solutions, Eq. (25) into the action, 

Eq. ( 24 ) , we arrive at the following contributions to the action: 



|did'^xvi/,,±(x,i) (^za, + ^ - t/|*,,±(x,t)p + ^^ *^,±(x,o + E a7tI^T(mF'* 

5a(x)<?;(x) 



E 



AA + C/|*3,±(x,t)P 



*<,,±(x,t)r<±(i)a^,±(t)- 



(26) 



In what follows we shall approximate Aa + [/|^g,±(x, t)p by A^, using the fact that the interatomic interaction is 
typically many orders of magnitude weaker than A^i. 



B. Eliminating the photon states 

Next, we integrate out the photon states, doing so in two steps: (i) by integrating out the environment modes to 
arrive at an effective action for the cavity photons, and (ii) by integrating out the cavity photon modes to arrive at 
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an effective action for ground-state atoms alone. To achieve step (i) we use the result of Caldeira and Castro-Neto 
[i.e., Eqs. (36) of Ref. j33j in the limit huJc ^ ksT] for the path integral over extracavity modes, and thus we identify 
the following contributions to the action that involve the cavity modes: 



Y^ / dtd'^x{fle"^^^gaix)n+ix,t)a*^^^{t) +h.c.) +Y^ / dt d'^x a^ _^(i)a^,+ (i)g„(x).g;;(x)n+(x, i) ~ {+ ^ ~) 



1 



a,P ' 



where n±{x,t) = '^j.{x,t)'i±{x,t) is the local atomic density. Note that we have dropped the atomic internal-state 
index g (i.e., we have made the relabeling ^g,± — >■ ^±). 

To achieve step (ii), we observe that the action is quadratic in the cavity photon modes, so they too can be integrated 
out, to produce the desired atom-only action which, for convenience, we express in terms of the classical-quantum 
(i.e., c — g) basis for the fields rather than the ± basis (see Sec. IV I: 



S°« = J dtd^X^l{x,t) (idt + ^ +m) *g(x,t) - t/*:(x,t)*;(x,i)[*,(x,t)2 + *,(x,t)2] +h.C. 



(28) 



-f ^Tr In M + / d'^x d'^x' dw dto' ^ 

a 

in which M is the matrix 

M(w,a;';a,^)= . ^^ 



and the quantities Di are defined as follows: 
1 



»^g^(x) 

Ai 



(ni(x,..) nMM)\^{u,,J-a,^)Y^ [l^^^^;'). 



n2(x',w') 



e/ '\x , ( Di(w'-i:j;a, /3) \ , , 



Di 



A, 



d'^xga (x) gp (x) n^ (x, w' - w) , 



(30) 



where the Keldysh components of the atomic density are 
given by 



ni(x,a;)= / dw' **(x,a;')*g(x, w - w') + h.c, (31) 
n2(x,a;) EE / do;' |**(x, w')\I'c(x, w - a;') 

+*;(x,c^')*«(x,^-^')}- (32) 

How one should proceed from here depends on the rela- 
tive magnitudes of k, Ap, and g'^N/Aj^. Our objective 
in this paper is to analyze the self-organization transi- 
tion in a multimode cavity. Physically, this transition 
is associated with the laser-cavity interference term, as 
discussed in Sec. |III[ and is most straightforward to an- 
alyze when El 3> Ec (i.e., when f2 3> 5a); moreover, as 
discussed in Ref. [S] (and as we shall show), the steady- 
state temperature of the system is proportional to Hk. 
Therefore, in order to explore self-organization at low- 
temperatures, it is natural to take Ac » K,g'^N/AA 
|34j . We may therefore expand InM and M~^ in powers 
of k/Ac and g'^N /{AaAc), thus arriving at a simplified 
atom-only effective action: 



r 



where the five terms — which are labeled by their corre- 
sponding coupling constants — respectively account for: 
5o, the kinetic energy of the atoms; 5*,;;, the A-periodic 
interaction caused by the scattering of photons between 
the laser and cavity modes; 5'^, the interaction due to 
the scattering of photons between cavity modes; Sjj the 
contact repulsion between the atoms; and S'k, dissipative 
processes due to the leakage of photons through the cav- 
ity mirrors. The terms have the following explicit forms: 



W2 

2A7 



Sn = j dojd'^x *: f cj + ^ + Ai I *g ^ 



h.c. 



(34) 



ni(x,w)n2(x',w) ni(x,w)n2(x',w) 



Ac +iK 



iC - IK 



S^ ^^AcY,] d^x d^x' S„ (x) S/3 (x) S„ (x') Ep (x') 

7ii(x,a;)n2(x',w) 



duj duj - 



SeS = Sq + Sq + S^ + Sjj + Sk, 



(33) 



{uj + Lu' — Ac -l- iK)(uj — Ac — ««:) ' 
Su = U f d'^xdt^*{xt)'^l{Kt)[i'^{xtf + ^q{xtf]+h.c., 

S,^ I dcu- r-- — - I duod'^xd'^x' ga{^)ga{-x!) 

(u) - Loc) + «: 7 

xi(^KAc coth(/iK//cBT) ni(xa;) ni(xcL'). 
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To streamline the notation we have introduced the cou- 
pling constants C = ^^^^^/(A^Ac) and ^ = ^''/(A^Ac). 
For the regime in which Ap is large compared with the 
other frequencies in the effective action (e.g., the typical 
atomic kinetic energy), a further simplification is possi- 
ble: the integrals over w and lo' can be expanded in a 
gradient expansion in terms of (1/Ac)9t. In what fol- 
lows, we shall keep only the zeroth order term in this ex- 
pansion, thus arriving at the following "instantaneous" 
forms of S(^ and S'^, in which we have expanded ni and 
n2 in terms of the atomic field operators: 



5^=C^ /'dirf''a;dVSc,(x)S„(x') (35) 

x*:(xt)^;(x'i)[*,(xt)*,(x't) + *,(xt)^,(x'i)]+h.c. 
5^=e 5] / '^^ ^""^ '^''^' S" W ^/9 W ^" (^') S/3 (x') 

X**(xt)*;(x'i)[*c(xt)*c(x't) + *,(xt)*,(x't)]+h.C. 

Note that if we neglect the effects due to S^, the nonequi- 
librium character of the theory would apparently disap- 
pear: the nonequilibrium laser- and cavity-mediated in- 
teractions 5*1;; and 5^ have precisely the same form in 
terms of Keldysh and time indices as the contact re- 
pulsion Sjj; thus, an effective equilibrium description of 
these terms should be possible. In the rotating-wave ap- 
proximation (see Sec. VI below and Ref. [7]), the laser, 
although a nonequilibrium element, only influences the 
density matrix by raising the system's "apparent" en- 
ergy eigenvalues by hiOL. As explained in the following 
section, we can use this fact to develop an effective equi- 
librium description of the self-organization transition, to 
which the effects of a nonzero 5k can later be added as 
a perturbative correction. 



VI. EFFECTIVE EQUILIBRIUM THEORY 

In this section we offer some considerations on self- 
organization in pumped cavities in general and, in par- 
ticular, on the quasi-equilibrium character of this tran- 
sition in the high-finesse limit. The fundamental differ- 
ence between the zero-temperature and nonequilibrium 
formalisms is the assumption — valid in the former case — 
that the initial state (in this case, the ground state) of the 
system without interactions evolves adiabatically into a 
pure energy eigenstate (in this case, the ground state) 
of the interacting system, provided the interactions are 
switched on sufficiently slowly. This assumption does 
not, in general, hold away from equilibrium, but there are 
certain nonequilibrium systems — e.g., an atom pumped 
by a far-off-resonant laser — for which it does hold: in 
this example, Fermi's golden rule implies that the laser 
does not stimulate real transitions between the atomic 
levels unless the laser is resonant with the atomic transi- 
tion. (For a broadened atomic level, Fermi's golden rule 



implies that real transitions are negligible as long as the 
laser-atom detuning, Ayi, exceeds the linewidth, 7, of the 
atomic transition.) 

In the present case, the excitation gap between the 
trivial ground state (i.e., the photon- free cavity and a 
uniform distribution of atoms) and the lowest excited 
state, which has photons and a A-periodic atomic density 
modulation (a "phonon"), consists of two parts: (i) the 
energy cost hl^c of adding a cavity photon and (ii) the 
energy gain due to the photon-phonon coupling, which in 
second-order perturbation theory would have the form 
{h'^{Vtg/l::^Af}/{fi'^Kl/2M}, where Kq = 27r/A. Thus 
the total gap must have the form: 



1 — const. X 



h^K^/2M 



(36) 



For the system at hand, the assumption of adiabatic- 
ity holds as long as the cavity is of sufficiently high fi- 
nesse, i.e., S > K. As we shall see in Sec. |XI[ the self- 
organization transition is weakly discontinuous, and the 
range of values of the control parameter Q over which the 
ordered and disordered phases coexist can be made larger 
than K. Under these conditions, the self-organization 
transition should be well-described by the effective equi- 
librium theory sketched in Ref. [T7] . Only for large values 
of K (i.e., bad cavities) would nonequilibrium effects have 
a chance of playing a leading role. For an alternative con- 
struction of such an equilibrium theory, which illustrates 
the relation between its quantum phase transitions and 
conventional equilibrium quantum phase transitions, see 
App.0 

The assumption of adiabatic switching implies that if 
interactions are turned on adiabatically in the distant 
past and turned off adiabatically in the distant future, 
an initial energy eigenstate would evolve into itself, up to 
a (physically unimportant) phase factor [35j . This con- 
sequence in turn implies that a single path integral can 
capture all the dynamical information; the second path 
integral in Eq. (13), which was demanded by the neces- 



sity of summing over all final states, would be superflu- 
ous because the final state would be known. In terms of 

Green functions, this implies that the components G^ 

and G \. become redundant; indeed, it can be proved di- 
rectly, in terms of Keldysh diagrammatics [3S] , that the 

expansion for G+-)- (and G ) is closed, containing all 

information about correlations and response. 

Formally, the case of nonzero temperatures is less 
straightforward because the adiabatic switching assump- 
tion is not available for mixed states, and therefore the 
derivation of the imaginary-time Matsubara formalism 
from the real-time Keldysh formalism is, in general, 
nontrivial. In this case, we are guided by the follow- 
ing consideration: if the fluctuation-dissipation theorem 
holds for the exact Green functions, the Matsubara for- 
malism is valid; in the present case, we can use the 
Keldysh diagrammatic technique (see Sec. VIII) to show 



that violations of the fluctuation-dissipation theorem are 
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small, and hence that the Matsubara technique is ap- 
proximately vahd. This is what one expects on physical 
grounds: the chief effect of the laser photons is to me- 
diate an effective atom-atom interaction, rather than to 
cause an energy flux. 

Our strategy in the next section, Sec. |VII[ will be to 
analyze the equilibrium critical behavior of the system 
at both T = and T 7^ 0, using established results from 
statistical mechanics and quantum field theory, and ne- 
glecting dissipative processes. After that, we shall return 
to the Keldysh formalism, in Sec. |VIII[ to reinstate the 
dissipative processes and explore their consequences. 



VII. QUASI-EQUILIBRIUM LANDAU- WILSON 
DESCRIPTION 

In this section we derive coarse-grained, Landau- 
Wilson forms of the atom-only action, Eq. (33), both 



for zero and nonzero temperatures, valid in the quasi- 
equilibrium regime. These Landau- Wilson actions are 
closely related to the one first introduced by Bra- 
zovskii |28| : we exploit this relationship in order to 
describe the impact of collective fluctuations on self- 
organization in multimode cavities. Finally, we discuss 
how the correlations of these fluctuations can be detected 
via the light emitted from the cavity. 

In an effective equilibrium theory, the prescription 
for going from the Keldysh to the Matsubara or zero- 
temperature formalisms is to trace, in reverse, the steps 
one would have taken to go from the equilbrium to the 
Keldysh formalism: i.e., keep the (-I--I-) component of the 
action and drop the Keldysh indices [31]. It is, further- 
more, convenient to reformulate the action in terms of an 
order parameter, i.e., a quantity that is zero in the uni- 
form phase and nonzero in the self-organized phase. The 
considerations of Sec. |III| suggest that the appropriate 
order parameter for detecting crystallization into cavity 
mode a should be given by 



Pait) = / (fxp{yi)ga{yi)- 



(37) 



The details of this procedure are dependent on the cavity 
geometry; we shall focus in this work on the case of the 
concentric cavity, as it is the most straightforward case. 



A. Mode structure of the concentric cavity 

A concentric cavity can be thought of as consisting 
of two mirrors that cover antipodal regions of the same 
sphere. Its mode structure is derived from the solutions 
of the Helmholtz equation inside the sphere. The effects 
of the edges of the cavity mirrors can be hard to com- 
pute accurately; the standard technique is to solve the 
Helmholtz problem approximately, by requiring the elec- 
tric field to vanish at the the mirrors edges. (As with 



the subsequent approximations that we shall make, this 
one works best for large cavities.) As the atomic distri- 
bution is quasi-two-dimensional, being confined near the 
equatorial plane, it breaks the spherical symmetry of the 
cavity, and it is therefore convenient to employ cylindrical 
coordinates (r, 6*, z). Thus, the planar dependence of the 
cavity modes takes the form Jm(^n^) cos(77i0), where Jm 
is a Bessel function, m is quantized by the requirement 
that cos(m6'o) = (see Fig.[T]), and n by the requirement 
that the field should vanish at the mirrors. The quanti- 
zation of solutions along the z direction (i.e., the pump 
laser axis; see Fig. fl]) yields a third mode index I. For 
large cavities, modes having a fixed value oil-\-m + n are 
frequency-degenerate. Because, as discussed in Sec. |ll) 
the free spectral range of the cavity is larger than the en- 
ergy scales relevant to self-organization, we can restrict 
ourselves to cavity modes having a certain fixed value of 
Aq = / -t- m -f n. 

The atomic density, which is quasi-two-dimensional, 
can be expanded over a similar set of mode functions, 
indexed by {m,n), provided one retains all such modes 
and not just those satisfying m + n = Aq. In addition, 
the boundary conditions on these mode functions are not 
in general the same as those on the cavity modes, as the 
atoms are confined by an external, confining laser field 
rather than by the cavity mirrors. For sufficiently large 
traps, however, this distinction is not expected to have 
important effects, and we shall neglect it. 

The cavity modes for which I = are expected to be 
favored for crystallization, as they have the highest am- 
plitude in the equatorial plane of the cavity, to which the 
atoms are confined. However, as we shall see, modes hav- 
ing I > are also of importance in setting the range of the 
effective atom-atom interaction, and thus in determining, 
e.g., the extent of the fluctuation-dominated regime, as 
well as the size and stability properties of droplets of the 
ordered phase [57] . 

In the (?Ti, n) basis for the atomic density, the order pa- 
rameter, Eq. (37), is given by p„m- Note, also, that pro- 



vided that the atomic density is spread out over a large 
number of optical wavelengths, the following asymptotic 
result holds: 



/ ^^ II, '^™'"' (^) ~ ^T.i "^^ ^JZi ' 



(38) 



In this expression, we have introduced the sign- 
insensitive Kronecker delta Sy^rrii = j 



E±'5e(-i)±'- 



Eq. (38 1 is exact for the angular (m) component, and 



holds approximately for the radial (n) component, if 
n > m; this is the regime of interest because modes for 
which m > n have large diffractive losses and do not 
couple to the atoms. Eq. (38 1 is closely analogous to 



momentum conservation, and simplifies the structure of 
quartic and higher-order terms in the action. 

It is useful, at this point, to specialize to two cases. The 
first addresses an ultracold gas of bosonic atoms, with- 
out contact interactions, which may or may not be Bose- 
Einstein condensed. We have dealt with this case, which 
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involves the introduction of an auxiliary field, in some 
detail in the Supplementary Information of Ref. [17j ; we 
revisit this case below in Sec. |VII B[ This approach can, 
in principle, be generalized to the case of an interact- 
ing gas (whether Bose- Einstein condensed or not), but 
is unwieldy for such systems, in which there are non- 
cavity-mediated interactions, as it involves the introduc- 
tion of multiple auxiliary fields. Thus, we treat the sec- 
ond case of interest, which is that of an interacting Bose- 
Einstein condensate (BEC) at temperatures well below 
the condensation temperature, using an alternative ap- 
proach that does not involve introducing auxiliary fields. 
Instead, we exploit the off-diagonal long-range order of 
the BEC and use a correspondingly modified form of the 
order parameter. In a BEC that is well below its conden- 
sation temperature, the atomic density factorizes to lead- 
ing order: (t/i'^ (r) ■0(i"')) ~ {'4'^i'^)) (^(i"'))- It is clear that, 
at temperatures much lower than the self-organization 



energy scale h'^KQ/2M, the low-energy modes — near the 
self-organization transition — are those corresponding to 
two widely separated regimes of "momentum" m + n, 
viz. m + n Ki and m + n Ki KqR/2'k. In the case where 
the system is Bose-condensed, one can use the presence 
of off-diagonal long-range order to exchange the order pa- 
rameter, Eq. (37), for the condensate amplitude (ipmn) — 



because self-organization then involves the macroscopic 
occupation of a mode with m + n — KqR/2tt [3S]. This 
exchange considerably simplifies the structure of the the- 
ory, and makes it possible to treat the effects of both the 
contact interaction between the atoms and the cavity- 
mediated interactions between them in a relatively trans- 
parent way. The final structure of the theory is, as we 
shall see, the same in both cases (i.e., the ultracold non- 
interacting system and the Bose-Einstein condensed in- 
teracting system). 



B. Ideal Bose gas 



In the ideal-gas case, we proceed as follows. We note that the cavity-mediated interaction term can be written as 
J dr^^^ QmnPmn{T) Pmn{T) , whcrc C,„„ is the cavity-mediatcd interaction favoring atomic modulation at wavenum- 
ber (m + n). As mentioned in Sec. VII A the coupling Qmn is to be considered as being peaked about modes obeying 
m + n = KoR/2tt. Next, we perform a Hubbard-Stratonovich transformation, which consists of introducing an 
additional, Gaussian functional integral into the partition function via the identity |32j 



exp ( / dT(„in Pmn{T) P-,nn{T) j = / D p„in CXp i-dr 



Pmnir) Pmnir) + 2— ^Pm„(r) Pmnir) ] , (39) 



in order to render the action quadratic in the ^ variables. The partition function can then be rewritten 



where 5" is given by 



S' 



dT / d''a;**(x,T) idr 



2M ^ h 



-2^-|-/5„„(r)S„„(x) *(x,r)-t- / dr^ 

mn / mn 



(40) 



f^ Sran 



PitnnXj' } PmnKj' ) • 



Provided the laser strength is below the self-organization threshold and the gas is Bose-condensed, the field operators 
^(x, t) can be expressed in terms of condensate and non-condensate parts as ^(x,r) = ^J Nq/A -\- $(x,t), where 
No{T) is the equilibrium condensate fraction at temperature T and A is the area occupied by the atoms. Transforming 
the Bose fields to the basis of (to, n) mode functions and Matsubara frequencies, in which the kinetic energy is diagonal, 
one has ^mn(wiy) — a//Vo '^m.o ^n,o (^ly^o + *&>Tm(w,y). Integrating out ($*, <&), one arrives at the action 



^"-iTrln(M) + ^ 



/ ^ y. Pmnu Pmn—u -^ '0 Pmnu l^-LVl )mnL',m'n'u' Pm'n'i/' 



(41) 



where the (infinite-dimensional) matrix M is defined by 



M 



mnv^m'n'v' 



-li^v 



h{m + nY 
2Mi?2 



^mn,m'n' 



2kBT sr^ . J- <■ 

■pqv" 



^v+v" ,}. 



(42) 
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Below threshold, so that self-organization is not present and (/5„ 
quadratic term in the action is then given by 



0, it is useful to expand M in powers of p\ the 



/ , PmnX P- 
mnX 






NokeT 



mn — 



lujx 






2fi2 



E 



1 






M 



h(m+n — {p+q))'^ 



P 



The last term can be usefully rearranged if one recalls that {ksT / h)Y^^^^^{iuJi, 2MB? m)^^ — N — Nq for a 

Bose- Einstein condensate [35]. (This statement also holds for a non-Bose-condensed gas, if one sets A^o — 0-) We now 
use the fact that the atoms in a Bose-Einstein condensate typically have energies that are small compared with the 
recoil energy to approximate the last term, and find that at low temperatures the quadratic part of action is then 
given by 



fceT 



E- 



' P-r 



1 



N 



Cn 



-ico^ + h{m + nY/2MW 



(44) 



in which the vanishing of the coefhcient for any {m,n^v) signals, at the mean-field level, the instability of the gas 
towards self-organization in the corresponding mode. (Note that self-organization is therefore possible only in the 
V = sector.) To enable an analytic treatment of the transition, we replace the (analytically inaccessible) exact form 
of C,mn by the convenient approximate form C,„„ = C,[l — x("t. + n — Ao)^], where the parameter 1/x represents the 
extent to which the coupling to I ^ Q modes is suppressed (see App. Ib] for further discussion of this point). This 
approximation captures the fact that the coupling of the atoms to the cavity modes is strongest for the modes that 
obey rn + n ~ Aq and otherwise simplifies the structure of the theory without making any drastic modifications to it. 
(For X = oo the I ^ modes are entirely suppressed; for x = the atoms couple equally strongly to all cavity modes, 
in which case there is no preferred lengthscale for self-organization.) 

Continuing with the expansion of ln(M) and M^^ in powers of /5, assembling the two contributions, and retaining 
the zeroth-order term in a gradient expansion, we arrive at the following form for the quartic-order term: 



{k^rfN 



rriin^, 



Prriini Pm^n^ Pm^ri^ Pm4n4 0^2 lUi ^^ Ui ■ 



(45) 



Finally, we make the rescaling p — > pyJhCjJiBTx] in terms of these rescaled fields, the action, to quartic order, assumes 
the following Landau- Wilson form: 



'LW 



E 



+ 



1 

_x 



1 



NC 



hKl/2M ~^ {hK^/2My 



+ {m + n - {KoR/27r)Y 



P'lnny P—ran—v 



(46) 



X^f9K^/4.^P 



/ , Pmrnivi P-m2ri2i'2 Pni^nsus PmiUiVi ''^ mi,0 ''j] t^ifl "ni+n2,n3+i' 



If r > 0, we can restrict ourselves to the v = Q sector of the order-parameter theory, as this is the only sector of the 
theory that plays an important role for thermal phase transitions. In this case, S'lw is an instance of Brazovskii's 
free energy [55]. We shall return to this "action" (which is, in effect, a free energy rescaled by ksT) and discuss its 
implications for the character of the self-organization transition in Sec. [VII D[ after first deriving a closely analogous 
effective action in the case of the BEC with contact interactions. The case of T = requires extending Brazovskii's 
analysis to quantum phase transitions; we discuss this case in Sec. |VII F| 

I 

three auxiliary fields. Instead, we exploit the fact that 
the "momentum" -space structure of the action simpli- 
fies considerably for the low-energy modes, which are the 
modes of interest because they are the ones that provide 
the critical fluctuations. This simplification is analogous 
to that which arises in Fermi-liquid theory owing to the 
constraint that all low-energy excitations must have mo- 
menta that are approximately equal to the Fermi mo- 
mentum |21| . In the present case, the operative con- 
straint is that all values of m -I- n must be either ap- 



C. Interacting BEC 



the interacting BEC, one begins with 
(33)], as in the previous subsection (Sec 



^< 



cff 



VIIBl, 



For 

[Eq. '_ 

but proceeds differently. It would be inconvenient to ap- 
ply the auxiliary-field technique to the present case be- 
cause the action contains three quartic terms, each hav- 
ing a different "momentum" -space structure; decoupling 
the action would therefore require the introduction of 
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proximately zero or approximately Aq = KqR/2t:. We 
denote the m + n w components of the atomic field by 
(/), and the to + n w Aq components by $. "Momentum 
conservation" i.e., Eq. (38), then implies that the follow- 
ing kinds of quartic terms are admissible: (i) four (p fields; 
(ii) two (f) and two $ fields; and (iii) four $ fields. (Pro- 
cesses involving three $ and one fields are suppressed 
because at least one of the <f> fields would have to have 
m > n, which would imply large diffractive losses |30).) 
For terms of type (ii) , it is clear that the only kinds of pro- 
cesses that survive to arbitrarily low energies are those 
in which the two $'s and the two 0's have the same val- 
ues of {m,n). [In principle, terms involving pairs ($,0) 
having the same value of n/m should also survive to ar- 
bitrarily low energies, but they can be shown to have 



negligible phase space, compared with the other terms 
mentioned.] Similarly, for terms of type (iii), the only 
sets of [m,n) that survive to arbitrarily low energies are 
forward- and backward-scattering processes, viz., those 
for which {mi,ni) — (7712,112) and (TO.3,n3) = (7714,714). 
The other processes are said to be "irrelevant at tree 
level" [2lj because they become progressively less impor- 
tant at lower energies. 

Applying the arguments just given to the three quartic 
terms in 5'eff, we find: (1) that at low energies the mode- 
mode scattering term S*^ is irrelevant for (/)'s (except for 
the term in which two of the incoming momenta are zero) , 
whereas it does survive for $'s; and (2) that the contact 
repulsion separates into three parts, and can be written 
as follows (note that all interactions are local in time; for 
this reason we have suppressed the time arguments): 



U 



/ J T^miTii r^m^n, 



sy^m2n2r'm4n4 



'HE' 



mriT^mn 



mirii 

Y. *m'„'*m'„' I + ( H '^mn'/'mn I I Y. *™'»'*™'"' I 
^m'n' / \mn / \m'n' / 



h.c. 



(47) 
(48) 



As in Eq. (38), we have used the sign-insensitive Kronecker delta, Ssnii = ^ ^_|_ S\2(-i)±mi- We can now integrate 
out the (j) modes, provided we first render the action quadratic in these modes; this can be done either by making 
the Bogoliubov approximation (see, e.g., Ref. [3S]) or, more generally, by exploiting the fact (which follows from 
Goldstone's theorem [32]) that the low- lying modes of a BEG are linearly dispersing phonons. If this is done, the 
action for the 6 fields assumes the form 



^0= / ^^im 



'-^ujmn T^—ujmn} 



inn rn'n' 



S T 



Tuj'in'n' 
^—uj'm'n' 



(49) 



This form of the action is, strictly speaking, only appropriate for zero-temperature; for T > the integral over to 
becomes the discrete sum fc^T ^^ . For compactness, we shall present only the expressions for T = 0, except when 
the two cases differ substantively. The blocks S and T are given as follows [note that both are diagonal in the {m, n) 
index; the appropriate delta functions have been omitted for compactness]: 



5„„(w,w') = 



him + nY 



SojLu' + Y, [U + C{in' + 112,71' +n)] Um'n'i^ - tj'), 



Tmn{^,^') = UtiqS^^^' + Y [C^ + C('77' + m, n' + 77)] ;/„/„/ (-w-w'). 



(50a) 
(50b) 



where no denotes the number-density of particles in the condensate, and 

We now integrate out the (p fields to arrive at the following effective action in terms of $: 

{Kr,M *™n(-^)) ( l^+ %# +W + V -W + V \ f $„„(^) 



(5Gc) 



S'[$,$* 



du^Y, 



-W + V 
TmW 



-W + V -tiO+^^^^^^+W + V \^*mni-^) 



r 



(51) 



where 



with A'o being the total number of particles in the con- 



Vmn — Uuq — CmnNo, 

yVmn = U{n -no) + I duj 



~CmnUno 
2 I 7- f ?»("!+"- A)^ 



(52a) 
,(52b) 
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densate. 

This action, Eq. (51), can, once again, be addressed 
by means of a Bogoliubov transformation. Upon per- 
forming such a transformation, we find that, for suffi- 
ciently low energies, the quadratic term S2 in the ac- 
tion can be expressed as follows in terms of the real field 
$(a;) = $(a;) + $*(-a;): 



S2^ldu;Y.^^ 



' hKlJ2MR^ ' 



(53) 
in which x again represents the coupling to modes having 
I > 0, but the control parameter r is now given by 



hKl 



2M V 2M 



^^Ku-CN 



(54) 



As for the term quartic in $, its leading-order gradient 
expansion is a contact interaction among the $'s, pro- 
portional to 



5'4=ilJ| Mwi ^ $A/,,Ar,,Ldj<5sj\/, few. '^(SWi). (55) 
i=l Mi,Ni 

The coefficient il accompanying this term receives contri- 
butions arising in two ways: (1) from terms in S^s that 
are quartic in $, and (2) via the integrating out of 4>. 
In principle, these each give rise to three terms, one as- 
sociated with each of the quartic terms the microscopic 
action, Eq. (33). However, owing to the "momentum" - 
space structure of this action, Sq contributes no terms of 
type (1), and S^ contributes no terms of type (2). On 
the other hand, Sjj contributes both types of term. Of 



these, the term of type (2) is proportional to the above- 
the-condensate density, which is expected to be small at 
the relevant temperatures; hence, this term is sublead- 
ing, compared with the term of type (1). Putting these 
facts together, we find that the quartic term, Eq. (55), 
has the coefficient 



il=U -^N 



u 



(56) 



At nonzero temperatures and near the self-organization 
transition, the sum over Matsubara frequencies is domi- 
nated by the a; = sector; under these circumstances, the 
action S'2 -I- 5*4 -I- • • • [Eq. ( |5l| )] has the same form as that 
for the noninteracting case, and they are both variants 
of the free-energy functional first analyzed by Brazovskii. 
At T = 0, the situation is somewhat different; we shall 
return to this case in Sec. IVII Fl 

That the actions have the same form is due to the 
phenomenon of "universality" near phase transitions: the 
structure of any theory sufficiently close to a critical point 
depends only on the symmetries of the order parameter 



and the free energy (or "action"). In the present case, 
there are two salient features: (1) that there is a strip 
of degenerate, low-lying atomic density modes around 
2iT[m -\- n) — KqR] and (2) that terms cubic in the or- 
der parameter, which would be allowed by symmetry, are 
forbidden because they would involve at least one mode 
having m > n; such modes have high diffractive losses 
and cannot, therefore, be effectively populated with pho- 
tons. These constraints are sufficient to force the action 
to have the above form near the phase transition. 



D. Classical Brazovskii transition 

In this section, we briefly outline our adaptation of 
Brazovskii's self-consistent analysis of the eponymous 
model |28] to the present setting. Brazovskii's analysis 
begins with the following free-energy functional in terms 
of the real order-parameter field ip{x) and its Fourier- 
space counterpart ipy^ = J dfix exp(ik • x) V'(x): 



F^ I d'-k i^k [n + (|k| - Kf] ^P-k+U f d'^x (V'(x))^ , 

(57) 
where the bare phenomenological parameters TZ and 14 
are, respectively, the control parameter for the transition 
and the interaction parameter. At the mean-field level, 
F is minimized for 7?. > by the uniform configuration 
'0(x) — 0, and for 7?, < by ^k having the nonzero 
value \J —TZlTlA for any one of the momenta k having 
magnitude /c^. 

In order to adapt Brazovskii's analysis to the present 
case, we must replace all instances of the momenta k 
by the sets of positive mode numbers (m, n) . In par- 
ticular, we must replace the expression (|k| — fcj,)^ by 
R^^{m + n — Aq)^, where R is the radius of the concen- 
tric cavity. Therefore, the low-lying excitations in the 
present geometry do not lie on a circular or spherical 
shell in momentum-space, as they do in the original Bra- 
zovskii problem; instead, they lie on a linear ribbon along 
m + n = Aq in mode-space (see Fig. [2]) . This difference 
does not affect Brazovskii's argument, except for some 
numerical factors of order unity, as we shall see below. 

The primary consequence of order-parameter fluctua- 
tions, to leading (i.e., one- loop) order, is to renormalize 
the bare parameters TZ and U in the free energy Eq. (57); 



thus, the fluctuation-corrected free-energy has the same 
form as Eq. (57), but with corrected parameters r and 
u instead of TZ and U respectively. For the former, one 
must evaluate the Feynman diagram in Fig. |3k; the result 
is that r is implicitly given by 



r = TZ + aKoUr- 



-1/2 



(58) 



in which the coefficient a in the "self-energy" term is a 
geometrical factor that, in the present case, is approxi- 
mately given by the expression 
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a) 



b) 



mnw mnoj mnco mnco mnco mnoj 
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FIG. 2. Dispersion relation for low-energy atomic excitations, 
i.e., those that approximately satisfy 27r(m + n) = KqR; as 
discussed in the text, the trough-like form of this dispersion 
enhances fluctuation effects. The inset shows a "top view" of 
the dispersion: the black line represents modes at the mini- 
mum of the trough, which exactly satisfy 27r(m + n) — KqR; 
self-organization results in the macroscopic occupation of one 
of these modes. 



a w 4\/2nmodcsA/L, 



(59) 



with riiiiodcs being the number of modes having finesse 
high enough to be populated. (As one might expect, fluc- 
tuation corrections are more important in cavities hav- 
ing a larger number of degenerate modes.) Note that, 
regardless of the sign of 7?., the corrected value of r is 
positive; hence the apparent second-order transition out 
of the disordered state is precluded by fluctuations. 

The leading corrections to lA are given, for generic val- 
ues of TO, n, to', n', by Fig. pp. For {m,n) = {m',n') one 
must also consider the diagrams in Fig. [Sfc. Summing 
up these series of diagrams, one finds that the corrected 
value u, for {m,n) = {m! ,n'), is of the form 



■ U 



1 - {U/r^/^) 
1 + (i^/r3/2) • 



(60) 



Evidently, u turns negative for sufficiently small r. 
Naively this would mean that the free energy becomes 
unbounded below; however, perturbative corrections en- 
coded in diagrams such as Fig. [Sji generate a positive 
six-point coupling, associated with a coupling denoted as 
w, which stabilizes the action and gives it the profile with 
multiple minima shown in Fig. [4] This profile for the free 
energy suggests that any phase transition that the system 
undergoes is likely to be first-order (i.e., discontinuous). 
To determine when the transition becomes energeti- 
cally favorable in equilibrium, one should compare the 
free energy of the uniform state with that of the possible 
self-organized states. (We shall revisit this point for the 
nonequilibrium case in Sec. VIII ) Let us consider, first, 



the state in which all the atoms are self-organized in a 
single mode, so that {ipmn) = A5mM ^n,n, the coefficient 



c) 




m n CO mnco 



mnco m n co 




d) 




FIG. 3. (a) Dyson equation for the self-energy at one loop 
order (i.e., the leading fluctuation correction to r). (b) A geo- 
metric series of corrections to the vertex (i.e., to u), which con- 
stitute the primary fluctuation corrections for (m' ,n' ,ij') ^ 
{m,n,Lj). (For the classical case, uj' = lj = 0.) (c) A geo- 
metric series of corrections to u that contribute only when 
{m' ,n' ,uj') « {m,n, ij). It is these contributions that change 
the sign of u, thus causing a first-order transition, (d) Higher- 
order vertices that emerge under coarse-graining. 



A being the amplitude of the order parameter. If $o is 
the bulk free energy of the disordered state and Fi that 
of the ordered state, one can formally write 



Fi — Fq 






dA 



SF d^„ 



mn 



dA 



(61) 



The motivation for this rewriting is that hmn = 5F/5ilirnn 
is the biasing field that would render a certain order- 
parameter configuration stable. Both the disordered and 
ordered states are locally stable at zero field; therefore, 
h should go to zero at both ends but should be nonzero 
between, so as to "drag" the system from one phase to 
the other. The advantage of integrating /i, as opposed to 
computing the free energies directly, is that one avoids 
having to compute terms in the free energy that are the 
same in both phases. The leading contributions to h can 
be written as follows: 



1 3 
6 



1 



■A^ 



1 



uAJ2{^^ 

mn 



Ymn) • 



mn Ymn 



(62) 



The last term should, in principle, be computed to the 
same order as fiuctuations have been computed in the 
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Order parameter 



FIG. 4. Schematic form of the free energy as a function of the 
order parameter, both above and below threshold, indicating 
how fluctuations change the character of the phase transition. 



any thermodynamic system with a continuous symmetry 
(and short-ranged interactions). This result is a con- 
sequence of the large phase space associated with long- 
wavelength fluctuations of the direction of ordering — in 
the present case, to fluctuations of the phase of the su- 
perfluid order parameter and/or the direction of order- 
ing |IX A[ Therefore, one would not expect an infinitely 
large sample to exhibit true long-range order at finite 
temperatures. This result is not, however, particularly 
relevant to the case at hand, for the following reason. At 
T = 0, the system is effectively three-dimensional rather 
than two-dimensional, owing to the additional dimension 
that corresponds to imaginary time; at sufficiently low 
temperatures, therefore, one expects the distances over 
which spatial fiuctuations destroy long-range order to ex- 
ceed the system size for relatively small systems such as a 
typical BEC [42]. Even at higher temperatures, one can 
suppress long-wavelength fiuctuations by using a system 
having multiple layers, so as to increase the effective stiff- 
ness against fluctuations of the order parameter. 



disordered state, viz. to one-loop order. At this or- 
der there are two diagrams that need to be computed: 
Fig. [3^ and Fig. [Sja. The corrections are substantially 
different for {m,n) = {fh,n) (i.e., the longitudinal com- 
ponent) and {m,n) ^ {m,n) (i.e., the transverse, or 
Brazovskii's anomalous, component). The longitudinal 
corrections have essentially the same form as those we 
computed in the disordered phase. The transverse cor- 
rections, however, diverge with system size, as a conse- 
quence of the Mermin- Wagner theorem (see the following 
section). Swift and Hohenberg [39] have shown that for a 
finite system in two spatial dimensions, these corrections 
are small as long as K^R < m~^/^, a condition that is 
met for sufficiently weak coupling [3D]. Neglecting these 
contributions, one finds that the free-energy difference 
between the disordered and ordered states is given by 



A$ 



^iV^^V^) 



1 
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where rA is defined implicitly via 



rA=n 
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and r is given via Eq. ( 58 ). Solving this pair of equations, 
one finds that A<i> changes sign — and the equilibrium 
phase transition therefore occurs — when TZ ~ —{aUY/'^. 



E. Relevance of the Mermin- Wagner Theorem 

A well-known result in the theory of phase transitions, 
the Mermin- Wagner theorem (see, e.g., Ref. jH]), states 
that long-range order is impossible in two dimensions for 



F. Quantum Brazovskii transition 

The quantum case of the Brazovskii transition, which 
occurs at T = 0, differs from the classical case in that 
the quadratic part 5*2 of the action governing it has the 
form given in Eq. (53), viz.. 
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(65) 
in which the frequency integration variable uj has been 
rescaled to absorb certain dimensionful factors. The 
quartic term in the action, Eq. ( 55 ) , also includes fre- 
quency integrals. The presence of these frequency inte- 
grals, absent from the classical case, changes the spec- 
trum of fluctuations. Qualitatively, this is because there 
are now two dimensions transverse to the critical sur- 
face; therefore, instead of a ribbon of critical modes, one 
must consider a tube. If one suitably adapts the Bra- 
zovskii diagrammatic procedure, one arrives at the fol- 
lowing implicit expressions for the fluctuation- corrected 
parameters: 
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In these equations, i? is a high-energy cutoff (which 
would be of order IS. a in the physical system) . It is tempt- 
ing to interpret them as follows: as 7?. — > — oo, we have 
that r — > 0; therefore, r ^ Beyi\){—\R\/aU), and thus 
r is always positive, although it does become exponen- 
tially small in the 7?, — >■ — oo limit. This would seem to 
suggest that criticality is not restored at zero tempera- 
ture. Furthermore, the fluctuation-corrected vertex has 
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the approximate form 



u : 
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(67) 



which suggests (cf. Sec. VII D) that metastabihty should 
set in when |7^| = aU\og{B/aU^). 

As these results are cutoff-dependent, however, one 
should investigate the quantum Brazovskii action us- 
ing a more systematic scheme than Brazovskii's, e.g., 
a renormalization-group scheme such as that devel- 
oped for the classical Brazovskii problem in Ref. [29] . 
We describe the appropriate quantum adaptation of 
this renormalization-group scheme in App. [Cj The 
corresponding renormalized values of the various (de- 
dimensionalized) parameters are shown in Fig. p^ as func- 
tions of the (de-dimensionalized) bare control parameter 
7?.. The fact that u first becomes negative for a smaller 
value of —TZ than the value at which r goes to zero indi- 
cates that the transition remains first order, and this is 
one of our main results. This result can, however, be de- 
duced on grounds that are more physically transparent, 
as we shall now discuss. 



Therefore, the renormalization of r in the present setting 
should be the same as that in the 0{p) model at large p. 
It is known that in (1 -I- 1) dimensions and for p > 2, the 
parameter r is always rendered positive by fluctuations; 
therefore, criticality is never achieved in the 0{p) model, 
and any phase transition that might occur in the 0{p) 
model — and hence the present model — must be first or- 
der. (This remark also applies to the classical Brazovskii 
case; the absence of criticality does not depend on the 
weak-coupling approximation that Brazovskii's analysis 
employed, but follows from the low-energy structure of 
the free energy.) 

Unlike the 0{p) model, which is isotropic in order- 
parameter space, the present model does undergo a tran- 
sition (which is first order), owing to the additional con- 
tribution to the quartic term involving all four 0's being 
equal. In order to show that such a transition is feasible 
free-energetically, one can turn to the renormalization- 
group scheme outlined in App. IC) 



H. Fluctuation-corrected threshold: summary of 
results 



G. Analogy with 0{p) vector model 

At low energies, the only four-point couplings that are 
relevant, and therefore survive under coarse graining, in- 
volve either forward- or back-scattering. Consider order- 
parameter modes {m,n) that satisfy 27r(m + n) = KqR; 
for these modes the quartic term takes the following form, 
in which time indices (which follow from locality in time) 
have been suppressed: 
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(68) 
It is convenient to introduce the notation: 6 = nim and 
rj= |Ao — (m-l-n)l, which together provide an alternative 
labeling of the mode (to, n). In terms of these labels, the 
action takes the form 
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If one had ignored the contribution from terms having 
four equal values of 9 (i.e., the Sg,g' term), this action 
would be an instance of the 0{p) model in (1 -|- 1) di- 
mensions, with p being the number of cavity modes sat- 
isfying the degeneracy condition 27T(m + n) = KqR. 
In the large p limit, terms having 6 — 9' do not con- 
tribute to the renormalization of TZ: the index 6' (or 
9) is summed over in the relevant diagrams, and the 
single value 9 = 9' contributes negligibly to this sum. 



We now list the fluctuation-corrected values of the 
threshold for self-organization, for the three cases dis- 
cussed in previous subsections, in terms of the physically 
relevant microscopic parameters. In each case, the quan- 
tity listed is the fractional change in the threshold pump 
laser strength, i.e., (fith — ^th)/^th'' ^'^^ ^^^ regimes 
in which our analysis is valid, this quantity is generally 
much smaller than unity. 

(i) For the ideal gas at T > 0, the corrected threshold 
for the Brazovskii transition is given by 
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This expression is somewhat simpler than that given in 
Ref. p. 7] for the absolute change in threshold, but is 
equivalent to it. (In the present expression, we have ex- 
plicitly included a factor a that is related to the number 
of cavity modes that couple appreciably to the atoms.) 
The fluctuation correction can be thought of as consist- 
ing of two components: (i) a is a geometric factor, intro- 
duced in Eq. ( [59| ; (ii) the other part of the equation is a 
measure of the number of particles per correlation area, 
in the sense that (as discussed further in Sec. |VIIIB| the 
quantity y'x/^o is analogous to a (far-from-criticality) 
correlation length. 

(ii) For the interacting BEC at T > 0, the shifted 
threshold is given by 
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FIG. 5. (a)-(c) Dependence of coarse-grained, fluctuation-corrected parameters on the bare control parameter 7?., which is 
related to the laser strength, for a fixed value of the bare parameter W. (The bars over the parameters signify that they have 
been rescaled as described in App. [c]) These results are obtained by integrating the renormalization-group equations derived 
in App.lC] Panel (a) shows the flow of the effective "control parameter" r (which remains positive). Panel (b) shows the flow 
of the effective interaction parameter u, which changes sign as discussed in the text. Panel (c) shows the flow of the emergent 
six-point coupling w. Finally, panel (d) plots the free energy as a function of the order parameter A for three values of 7?., viz. 
—5.3 (thin solid line), —5.35 (dashed line), and —5.4 (thick line). The first-order phase transition takes place at 7?. « —5.36. 
These results are interpreted in terms of microscopic parameters in Sec. |XI[ 



Note that this shift explicitly depends on temperature. 
At sufficiently low temperatures this shift picks up fur- 
ther corrections, and eventually crosses over to the quan- 
tum Brazovskii result discussed in Sec. |VII F[ which, in 
terms of microscopic parameters becomes 
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In the fourth case (viz., that of an ideal gas at T = 0) 
power-counting suggests that the transition should re- 
main second order, even after cavity-mediated interaction 
effects are included. We have not discussed this regime 
in detail in the present paper; we plan to address its 
properties and achievability in future work. 



I. Signatures of criticality 

The Brazovskii transitions, both classical and quan- 
tum, being first order, do not exhibit the power-law de- 
pendencies (e.g., of the fluctuation correlation length. 



the order parameter, and its susceptibility) commonly 
associated with continuous phase transitions. However, 
as the Brazovskii transitions are only weakly first-order 
(i.e., they involve small discontinuities in the order pa- 
rameter), the influence of fluctuations on the atomic and 
optical correlations should be experimentally accessible. 
In particular, the fluctuation corrections to the density- 
density correlation function — i.e., Eq. [58] — should be ex- 
perimentally detectable in one of two ways. A straight- 
forward way to detect these correlations is to release 
the atoms from the trap, and analyze the correlations 
in the noise of the spatial density profile of the atomic 
system. These correlations can then be related to the 
density-density correlation functions of interest by means 
of the scheme described in Ref. 43J, which involves the 
post-processing of absorption images. An alternative 
approach — unique to the cavity QED setting, and prefer- 
able in that it does not automatically destroy the BEC — 
is through the correlations of the light emitted from the 
cavity. At weak coupling, the intracavity-^Aioiow corre- 
lations are directly related to the atomic density correla- 
tions, as follows. The full action (up to Gaussian order) 
is of the form loo^ a + {a) + a)p + pGp, where G is the 
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atomic correlation function given by Eq. 58 one can in 
principle integrate out p, thus arriving at the relation 
(a^a) ^ l/{uj ~ G^^). The fluctuation corrections to G 
are therefore manifest in the correlations of the emitted 
light. 

The weak-coupling approach, just given, has a serious 
limitation when it comes to describing quantum fluctu- 
ations: as discussed in the following section, hn acts, in 
some ways, as an effective temperature for the atoms in 
the cavity. Therefore, quantum effects are typically cut 
off by decoherence on a timescale comparable to 1/k. 
However, 1/k is also the timescale on which photons leak 
out of the cavity. It would seem to follow, therefore, 
that effects associated with coherent, quantum fluctua- 
tions take place on timescales too rapid to be detected 
via the leakage of light through the mirrors. This analy- 
sis, however, neglects the existence of cavity modes that 
couple relatively weakly to the atoms (so that they barely 
affect the effective temperature) and have considerably 
lower finesse. It is plausible, then, that the correlations 
of these modes can be used to probe the dynamics of the 
quantum fluctuations of the remaining degrees of freedom 
of the system. 



VIII. NONEQUILIBRIUM EFFECTS AT THE 
BRAZOVSKII TRANSITION 



In this section we reinstate the dissipative effects due 
to the cavity photon leakage rate k, and consider their 
impact on the self-organization transition. As we shall 
discuss, the departure from equilibrium implied by a 
nonzero value of k has three kinds of consequences: 
(i) it cuts off critical fluctuations, (ii) it affects the 
timescale on which the system is able to escape from 
a metastable state, and (iii) it modifies the dispersion of 
long-wavelength excitations in the ordered state. In this 
section we consider effects (i) and (ii); our discussion of 
effect (iii) is postponed to the next section, in which we 
discuss the properties of the ordered state. 

To accomplish this reinstatement, we follow the stan- 
dard prescription for expressing an effective equilibrium 
action in the nonequilibrium formalism (see, e.g., Sec. 4.7 
of Ref. [31]) and then augment the theory with the term 
5^ in Eq. ( |33[ ). Thus we arrive at a theory containing a 
copy of the equilibrium Brazovskii action involving fields 
on the ± contours, coupled to one another via Sk,, the 
consequences of which we shall now address using per- 
turbation theory. 



A. Critical effects 

In this section we focus on the case of the interacting 
BEC at zero temperature (cf. Sec. VII F), as the analysis 
is most transparent for this case. We begin with the 



dissipative term. 



S. 



duj 



d X d x' ga(x) ga{x.') 



{lu - Ac)2 + ^2 
x^CkAc coth{hK/kBT) ni{xuj) ni(xaj), (71) 

and re-express the atomic field 4* in terms of the conden- 
sate and non-condensed parts, as follows: 
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(72) 
(73) 



In terms of this decomposition, the primary quadratic 
contribution to the complete action arising from S^ then 
becomes 
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The most salient feature of this term is that its prefac- 
tor is frequency- independent. For a system to be at a 
quantum critical point, it is necessary that the prefactor 
vanish as w — > (see, e.g., Ref. [44| and our App. [d| . 
When — as in the present case — this condition fails, the 
collective dynamics on sufficiently long timescales is clas- 
sical. Indeed, the term Eq. (74 1 is formally identical to 
the term 
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(75) 



that arises for an otherwise isolated complex field $ cou- 
pled to an equilibrium thermal environment that is at 
temperature T [see, e.g., Eq. (66) of Ref. [31]]; one can 
therefore regard the coefficient Hk/kB = HkCN / {Ac^b) 
in Eq. (74) as an effective system temperature. In par- 



ticular, quantum correlations on timescales longer than 
1/k are washed out by the decoherence arising via the 
leakage of photons from the cavity: this effect is analo- 
gous to the decoherence due to a finite temperature that 
is known to occur near to a quantum critical point. 

This point is — in principle — immaterial, as interaction 
effects preclude criticality regardless of the value of k 
(as a consequence of Brazovskii's argument); however, 
at T = the fluctuation-corrected equilibrium control 
parameter is exponentially small, behaving as exp ( — 
\R\laU), and therefore it should be possible to tune the 
system close enough to criticality that the nonequilibrium 
suppression of criticality due to dissipation is observable. 

We now turn to the issue of influence of the nonequi- 



librium terms given in Eq. ( 74 ) on the effective quartic 
interaction vertex. At tree-level, the interaction vertex 
only couples terms having Keldysh indices cqqq or cccq. 
However, the nonequilibrium terms generate an effective 
ccqq vertex, via the Feynman diagrams shown in Fig. [6| 
This vertex brings the factor 
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FIG. 6. Contributions to the nonequilibrium vertex having 
external indices ccqq. For an introduction to the Keldysh 
diagrammatic notation see Ref. |31) . 



in which u is the fluctuation-corrected equilibrium vertex 
(see Sec. VII D). The most notable feature of this vertex 
is that it does not change sign when u does. In the regime 
considered in the present work, we have that k <C Ac, 
and hence u provides a subdominant correction to u, and 
therefore cannot prevent the net vertex from changing 
sign, signaling a first-order transition. It is possible, how- 
ever, that in the opposite regime, in which n 3> Ac, this 
correction term would dominate; in this regime, this term 
might be capable of preventing the Brazovskii transition 
from taking place at all. 



B. Nucleation and state selection 

In this section we address the dynamics of the emer- 
gence of self-organization associated with the Brazovskii 
transitions, classical and quantal. As these transitions 
are first order, one expects them to exhibit regions of 
two-phase coexistence, in which some parts of the cloud 
have self-organized and others have not. The time inter- 
val that the system takes to approach the steady state, 
in which the entire system is self-organized, depends on 
the energetics of critical droplets of the minority phase, 
which determines their nucleation rate. There are three 
regimes of interest, distinguished by the primary mecha- 
nism responsible for fluctuations: (i) near zero tempera- 
ture in an isolated system (i.e., in the regime where both 
the system temperature and h are smaller than [/), quan- 
tum tunneling is the primary cause of nucleation; (ii) at 
temperatures that are high compared with U , in an essen- 
tially isolated system (ksT ^ k), it is thermal activation 
that is the primary cause; and (iii) near zero temperature 
in the far-from-equilibrium regime (i.e., when k exceeds 
the temperature and U) , nucleation is primarily triggered 
by extrinsic force noise that originates with fluctuations 
in the photon population in the cavity. Owing to the for- 
mal analogy between k and T, discussed in Sec. |VIII A| 
and App.pj cases (ii) and (iii) can be treated by similar 
means. 

In all three cases, an essential ingredient is the energy 
barrier for thermal (or quantal) nucleation. In many set- 
tings this would be easy to read off from the tree-level 
(i.e., mean-field) Landau free energy (or action); how- 
ever, in the case at hand, the Landau free energy does 
not predict a first-order transition at tree level. On 
the other hand, it is not prima facie legitimate to use 
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the fluctuation- corrected free energy that was calculated 
above, as this incorporates fluctuations on all length 
scales, including those larger than the droplet itself. In 
general, therefore, one must follow a procedure like that 
due to Hohenberg and Swift [21], in which only fluctua- 
tions on length-scales smaller than the droplet diameter 
are self-consistently integrated out. 

Motivated by the wish to obtain analytical results, we 
focus, in the present work, on nucleation kinetics in the 
regime in which the corrected free energies of the self- 
organized and uniform phases are sufficiently similar that 
the critical droplet size is comparable to the system size. 
In this regime, the error incurred by using the fully renor- 
malized bulk free energy is expected to be relatively small 
(as, in this case, most of the renormalizations should al- 
ready have taken place) , and one can legitimately use the 
fully renormalized parameters computed in Sees 
and lynFl 



1. Classical nucleation 

In the classical case, one needs to compute the free- 
energy barrier to the nucleation of a critical droplet, 
i.e., the smallest possible droplet intrinsically capable of 
growing until it encompasses the entire cloud. The stan- 
dard procedure for doing this is to identify appropriate 
saddle-point configurations of the effective free energy. 
As these saddle-point configurations involve droplets of 
the ordered phase immersed in the uniform phase, one 
is interested in spatially varying configurations *&(x), 
and the position coordinates are therefore the appropri- 
ate ones to consider. In these coordinates, the Euler- 
Lagrange equation corresponding to the Brazovskii free 
energy takes the approximate form: 

e^(V2+x2)'*W + [r-2|u||$(x)p-|-3u;|$(x)|4]$(x) = 0, 

(77) 
in which ^g ^ ^/x/^o is an effective healing length for 
the crystalline order parameter. In order to find the 
saddle-point configuration, we look for solutions of the 
Euler-Lagrange equation that obey the boundary con- 
ditions that $ = near the edge of the cloud and 
^(x) ~ A'E.mn{x) near the middle of the cloud, where 

A = ±J {2u + V4u^ — 12rw)/(!>w is the value of the or- 
der parameter in the ordered state, and (m, n) is the 
mode into which the atoms are self-organized. In the 
conventional Brazovskii problem, the difference in en- 
ergy between a configuration including a droplet of area 
A and the uniform state is given by 
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in which Fq {— 0) and Fa are, respectively, the free en- 
ergy densities of the uniform and self-organized states; 
CTj^ (cTji) is the energy cost of an interface perpendicu- 
lar (parallel) to the lamellae; and L±^ (Lu) is the length 
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of the interface that lies perpendicular (parallel) to the 
lamellae. The concentric-cavity geometry differs from 
that considered in Ref. [53] in that it is not translation- 
ally invariant in its radial direction: both the atomic den- 
sity and the mode functions depend on x and so, there- 
fore, do the parameters in the Brazovskii model. As our 
purpose in the present section is to focus on order-of- 
magnitude estimates, we shall neglect this complication. 
A further difference between the present case and the 
conventional Brazovskii problem is that the modes are 
checkerboard-shaped rather than lamellar (see Ref. [T7], 
and also Sec. IX A ) . Thus, a generic interface has some as- 



pects of both transverse and longitudinal character, and 
the optimal droplet shape varies from mode to mode. In 
what follows, we focus on the (physically most relevant) 
modes for which tti <C n, and consider droplets of the 
form sketched in Fig. [7 

Returning to Eq. ( 78 ) , we see that the free-energy dif- 
ference per unit area, AF is of order 
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in which re is the value of the control parameter at 
which the equilibrium transition occurs. As argued in 
Ref. 35] , the width of a longitudinal interface is approx- 
imately ^o(^oCo/2r), whereas that of a transverse inter- 
face is ^0- Near the phase transition, r « {aU)"^^^, which 
is small relative to i^o^o for weak coupling U; hence, 
longitudinal interfaces are larger than transverse ones. 
Furthermore, the interface energy cost per unit area is 
given by the interface width multiplied by the quantity 
/o = ■u'^/(27w^), which is related to the curvature of the 
free-energy landscape about the minimum corresponding 
to the self-organized state. Thus, the total interface en- 
ergy cost is given by the expression 
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Because transverse interfaces cost less energy at weak 
coupling, the optimal droplet shape is needle-like, as 
shown in Fig. [7j 

To determine the free-energy barrier for droplet nu- 
cleation, we use the well-known Wulff construction (see, 
e.g., Ref. |45]); in the present case, this amounts to con- 
sidering bubbles of dimension oc an in the perpendicular 
direction and oc a± in the parallel direction. Thus, the 
free energy of a configuration including a droplet is: 



-Fdrop — ^F£,o 



2-f^oCo 



2r 



(82) 



The nucleation rate follows directly, being given by 
r exp(— i^drop/fcsT') or, in the case that nucleation is due 
to external noise [i.e., case (iii)], by rexp(— i^jrop/^^)- 
For a discussion of the relevant experimental parameters, 
see Sec. IXH 



2. Qualitative features of the defect morphology 

Near coexistence (i.e., for r ^ re), the critical droplet 
is arbitrarily large and — even for a highly anisotropic 
droplet — the energetic cost of the longitudinal interface 
(which scales linearly with bubble size) becomes greater 
than that of introducing localized defects, arranged so 
that the surface of the bubble is made as transverse as 
can be. A possible arrangement of such defects is shown 
in Fig. [7] It was argued in Ref. ^3] that such defects 
are energetically favorable only for \{TZ — TZc)/TZc)\ < 
([//^A^)^*^/^^, which is a narrow range compared with the 
thermal-fluctuation-dominated regime, which obtains for 

|(7^-7^,)/7^e)|<l. 



3. Quantum tunneling 

In the case of quantum tunneling rather than ther- 
mal barrier crossing, the argument of Sec. |VIIIB l| must 
be modified in two ways. First, the expression for the 
tunneling rate should be given by the form wq exp(— 6*0), 
where ujq is a characteristic collective frequency (e.g., be- 
ing proportional to the value of the renormalized control 
parameter r in the disordered phase) and 5*0 is the ap- 
propriate instanton action (see, e.g., Ref. [32] )■ A crude 
approximation to Sq is the product of the width (which 
is of order A) and height of the energy barrier; thus the 
difference between the initial and final values of the order 
parameter would act as an effective inverse temperature. 
Second, one must use the coarse-grained values of r, u, 
and w from the quantum rather than the classical model, 
i.e., from Sec. [VlTFJ rather than Sec. [VITD| 
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where 7 is a variational parameter that sets the overall 
scale of the droplet. By minimizing i^drop with respect 
to 7, we see that the critical bubble is that for which 
7 = TcHTc ^ ^); thus, the free-energy barrier for the 
thermal nucleation of droplets of the ordered state is 



IX. PROPERTIES OF THE CRYSTALLINE 
STATE 



In this section we examine various properties of the 
crystalline state exhibited by the coupled atom-light sys- 
tem, including its basic structure, elementary excitations, 
and topological defects. 
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FIG. 7. (a) Wulff droplets, corresponding to the TEMoo 
mode (j), and to a higher-order mode (ii), respectively. The 
droplets should become less anisotropic (i.e., less "needle- 
like") for higher-order modes; it is, however, possible that 
the optimal droplets in these cases have more complicated 
shapes, (b) Defected droplets, which are favored for r « r^, 
as discussed in the text. For these, the energetic cost of intro- 
ducing defects inside the droplet is outweighed by the increase 
in the fraction of the interface that is transverse. 



tical mode structure, and therefore on the possibilities 
for atomic crystallization. Instead, in the concentric cav- 
ity, the ordered states follow the two-dimensional opti- 
cal mode patterns, which may be visualized as distorted 
checkerboard patterns, as shown in Fig. [8] Locally, the 
atomic density selects amongst the cavity modes by crys- 
tallizing into the "odd or even squares" of the selected 
mode (cf. Fig. [s] and the discussion in Sec. III). In phys- 
ical realizations, the corresponding states of crystalliza- 
tion would not be exactly degenerate, as optical modes 
having stronger angular variation (i.e., larger m) are of 
lower finesse; on the other hand, repulsive interactions 
have a stronger impact on atoms that are crystallizing 
into modes of lower m 46 . Such effects can readily 
be accounted for within our model, via the introduction 
of fictitious fields that would bias the system towards 
crystallizing into certain modes. In practice, the most 
straightforward way to include such effects is by adding 
the terms 
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FIG. 8. Elementary excitations of the self-organized state in 
the concentric cavity, (a) Domains that have self-organized 
into distinct modes can be separated by analogs of grain 
boundaries (left half of panel) or by continuous textures (right 
half of panel), (b) Excitations that are anal ogous to the splay 
mode in smectic-A liquid crystals (see Sec. IX B I. Lines indi- 
cate nodes of the cavity electromagnetic field. The curved 
wavefronts along the radial direction have been drawn as flat 
lines to emphasize that the sketched feature is small-scale, 
relative to the size of the cavity. 



A. Basic properties 

1. Concentric cavity 

Soft condensed matter systems that undergo the Bra- 
zovskii transition commonly exhibit one-dimensional, 
lamellar patterns. The present realization does not, ow- 
ing to the influence of boundary conditions on the op- 



to the Brazovskii action, which would raise (or lower) the 
threshold laser power in a mode-dependent way. 



2. Other multimode cavities 

A well-known example of a multimode cavity is the 
confocal cavity, in which all the even TEM modes are 
degenerate i3_Q,. (It is also possible to make multimode 
cavities in which every p*^ TEM mode is degenerate.) 
These cavities have the practical advantage over the con- 
centric cavity that their stability criteria are easier to 
fulfill (e.g., their mode structures are more robust with 
respect to mirror misalignment). To the extent that it is 
legitimate to think of such cavities as each having a con- 
tinuous family of degenerate modes (i.e., provided they 
possess a large number of modes that are both degen- 
erate and not heavily suppressed by diffractive losses), 
the self-organization transition in these cavities should 
belong to Brazovskii's universality class, and our analy- 
sis of the transition itself should extend to these models. 
Where confocal cavities are likely to differ from concen- 
tric ones is in the geometry of the ordered states and of 
their defects, which is much more involved in the confocal 
case because of its less-evident symmetry structure. 

A feature common to most multimode geometries is 
that it is possible to tune the system across the point at 
which the modes are degenerate by gradually changing 
the mirror spacing. Thus, one can explore the crossover 
between the multimode physics discussed in the present 
work and the single-mode physics realized in Ref. ^16) . 
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FIG. 9. Case of the large-solid-angle concentric cavity, in 
which the atom-light system possesses a continuous symmetry 
associated with the relative phase between the +m and —m 
components of each mode function. This symmetry, when 
broken by the self-organized atomic cloud, leads to the ex- 
istence of both phonon excitations (shown in the left panel 
of the figure) and true edge dislocations (shown in the right 
panel of the figure). 



B. Phonons and nonequilibrium elasticity 

Manifestly, the concentric cavity geometry does not 
possess translational invariance in the radial direction; 
hence, there are no translational Goldstone modes in the 
radial direction. The geometry that we have primarily 
considered in this paper does not possess translational 
invariance in the angular direction either, because of the 
hard-wall boundary conditions that we imposed on the 
mode functions at the edge of the cavity, when comput- 
ing the mode structure (see Sec. VII A). For the usual 



experimental situation, in which the cavity mirrors cover 
a relatively small solid angle, this is the relevant case. 
In the opposite regime, in which each cavity mirror oc- 
cupies most of a hemisphere, one would essentially re- 
cover translational invariance in the angular direction; 
concomitantly, there would be phonons, corresponding 
to the "rippled" atomic arrangement shown in Fig. |9] 
For large cavities, such excitations would have a linear, 
phononlike spectrum, with a speed of sound related to 
the order parameter for crystallinity. 

The concentric cavity geometry does, however, pos- 
sess an analog of rotational invariance, in that the en- 
ergy is unchanged if one reorganizes the crystallization 
of atoms in mode {m, n) into crystallization in a degen- 
erate mode {m',n'). For a large cavity, having many 
modes, this symmetry is effectively continuous. Conse- 
quently, there are low-energy excitations involving the 
gradual variation of m and n across the cavity; an ex- 
ample of such an excitation is sketched in Fig. [8| The 
physics of these layer- wandering excitations is analogous 
to that of the splay mode of smectic liquid crystals |1T] ; 
as in the liquid-crystal case, the effective elastic energy 
for these excitations takes the Landau-Peierls form: 



where 9 (= n/m) parametrizes the macroscopically oc- 
cupied mode, and the ellipses indicate terms involving 
higher powers of the gradient operator. Near the tran- 
sition, the wandering rigidity Ki is proportional to the 
square of the equilibrium order parameter and also to the 
fourth power of the healing length ^o hi Eq. (77). 



We note in passing that the nonequilibrium character 
of the phase transition affects the spectrum of phonon 
and wandering excitations at very long wavelengths. 
For example, the effective retarded Green function for 
phonons, which in general has the form 



GR{uj,m,n) 



up- — Kcff{m + n)^ + 2iu)k 
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has purely imaginary poles when y/K^g {m + n) < k; 
modes satisfying this criterion are diffusive rather than 
propagating. This idea, which was discussed in Refs. [13 
|48] in the context of excitonic condensates, has nontriv- 
ial consequences for, e.g., the spatio-temporal decay of 
correlations in sufficiently large systems. 



F,i = Ki I d?x{y\ef 



(84) 



Defects 



In addition to the low-energy splay excitations (which 
are analogous, in some ways, to phonons), the ordered 
state can also have gapped excitations or defects, which 
in the present case are analogous to grain boundaries (see 
left panel of Fig. 8]) . When the ordered states on the two 
sides of such a boundary are of opposite parity, as in 
the figure, the boundary wipes out a fraction of a row 
of crystalline order. Its energetic cost is therefore ap- 
proximately L(7±^ , where L is the length of the boundary 
and (Tj_ is the interface energy discussed in Sec. |VIIIB| 
These defects are analogous to conventional topological 
defects in the sense that, for certain configurations of the 
order parameter at the boundary of some region in the 
cavity (e.g., at the edges of the cavity), the system is 
forced to have at least one defect somewhere inside this 
region. These defects are not, however, directly related 
to the existence of continuously broken symmetries in the 
system. 

Other kinds of topological defects might also be realiz- 
able. For instance, in the large-solid-angle case discussed 
above in Sec. IX B genuine edge dislocations, of the kind 
illustrated in Fig. 9] may arise. Another possibility is a 
texture of the kind sketched in the right-hand panel of 
Fig. [8} such a texture would be analogous to a closed 
lamella in the conventional Brazovskii case. It is not 
clear, however, that such textures are experimentally fea- 
sible^ as they would require that the system self-organize 
into a high-TO mode in some region of the cavity. 
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FIG. 10. Proposed scheme for detecting supersolid order, 
(a) Profiles of two cavity modes: Mode 1 (into which the 
atoms self-organize) and Mode 2 (which can be used to de- 
tect phase coherence, as discussed in SecH). The two modes 
are degenerate; Mode 2 possesses more nodes along the z di- 
rection (i.e., perpendicular to the plane of the figure). The ± 
signs describe the phases of the electromagnetic fields in the 
two modes relative to some reference (e.g., the pump laser) 
in various regions of the cavity, (b) Atomic configuration in 
which the atoms emit constructively into Mode 1 and destruc- 
tively into Mode 2. In the insulating phase, this is the typical 
configuration, as the number of atoms per site is fixed; hence, 
there is suppressed emission into Mode 2. (c) Atomic con- 
figuration in which the atoms emit constructively into both 
Mode 1 and Mode 2. Such configurations, which involve mul- 
tiple occupancy, occur in the superfluid phase but are sup- 
pressed in the insulating phase; hence, the amount of light 
emitted into Mode 2 is a measure of superfluidity. 



SUPERSOLID ASPECTS OF THE 
SELF-ORGANIZED STATE 



The spatially ordered state of a BEC in a multimode 
cavity is a "supersolid" in the following sense. It pos- 
sesses emergent forms of both crystalline and super- 
fluid order: i.e., it spontaneously breaks two continu- 
ous symmetries, the rotational invariance of space (to 
the extent that the cavity admits an effectively contin- 
uous family of modes) and the UiX) invariance asso- 
ciated with the phase of the condensate wavefunction. 
The properties, and even the existence, of supersolids 
have recently been active issues in condensed-matter re- 
search ini [m HSl HSl HH] . Amongst traditional condensed- 
matter systems, the primary candidate for exhibiting su- 
persolidity is solid *He, which was conjectured to have a 
supersolid phase in the late 1960s iSO, 51. Shortly there- 
after, Leggett [55] predicted that a supersolid would ex- 
hibit "nonclassical rotational inertia" when rotated suffi- 
ciently slowly, owing to the quantization of angular mo- 
mentum of a rotating superfluid. Evidence for this phe- 
nomenon was reported in Ref. fl)'; the interpretation of 



this and subsequent experiments is, however, still contro- 
versial. It has been proposed, for instance, that rather 
than indicating bulk supersolidity, the missing moment 
of inertia arises because of superfluidity in dislocation 
cores [121 [23] , because of elastic effects arising from the 
presence of ^He impurities [21], as a by-product of glassi- 
ness "^j and so on. 

One must distinguish between two types of questions 
regarding supersolids: (i) whether solid ^He, or any other 
neutral substance having "realistic" interparticle inter- 
actions, is supersolid in the sense defined above, and 
(ii) what properties a supersolid would possess, should 
one exist. In the condensed-matter context, attempts 
to address question (ii) have been vitiated by the un- 
certainty about whether the substance being studied is 
in fact supersolid, whereas attempts to address question 
(i) have been hampered by imperfect understanding of 
the characteristic phenomenology of supersolids. The ad- 
vantage of ultracold atomic realizations of supersolidity, 
such as the present one, is that one can explore ques- 
tion (ii) without first addressing question (i), as the ex- 
istence of both superfiuid and crystalline order is rela- 
tively easy to establish. The existence of crystallinity 
can be deduced via the superradiant emission of light 
(see Sec. VIII), whereas that of superfluidity may be ex- 
plored using a range of standard techniques (see, e.g., 
Ref. [ID). 

Furthermore, it is in principle possible in the present 
setting to test for supersolidity in a manner that enables 
one to distinguish between scenarios in which the bulk of 
the sample is a supersolid, and those that involve phase 
separation of some kind — e.g., scenarios in which BEC 
is restricted to, e.g., dislocation lines. This can be ac- 
complished via an analysis of the spatial correlations of 
the light emitted from the cavity. Such a technique is 
an adaptation of that developed by Ref. [53] to explore 
the superfluid-insulator transition in an optical lattice. It 
uses the fact that the number of particles per site is not 
fixed in a superfluid; therefore, even if the emission into 
a particular cavity mode is zero on average because of 
destructive interference between the contributions from 
even and odd sites, local atomic number fluctuations 
would render this destructive interference imperfect and 
would lead to a nonzero photon population, which can 
be detected in the light leaking out of the cavity. This 
idea is sketched in Fig. [TOJ for calculational details, we 
refer to Ref. 1551. 



A. Coupling the superfluid order parameter to the 
solid order parameter 

A well-known manifestation of supersolidity is the non- 
classical behavior of the the moment of inertia [S2] , which 
results from the requirement that the macroscopic wave- 
function be single valued. This effect can be explored 
directly in the present setting, e.g., by imparting angular 
momentum to the BEC via an auxiliary laser beam that 
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carries orbital angular momentum. Furthermore, one can 
study the implications of the presence of crystalline order 
for the superfluid transition [SI] by increasing the den- 
sity of the atomic cloud in the cavity until it undergoes 
Bose-Einstein condensation. A third possibility is to dis- 
place the superfluid from the center of the overall dipole 
trap and observe its relaxation [55|; this should provide 
information about the coupling between the Anderson- 
Bogoliubov mode of the superfluid and the excitations 
(both phononlike and topological) of the crystal. In par- 
ticular, it should be possible to tune the cavity geometry 
across a multimode geometry by adjusting the mirror 
spacing, as discussed in Sec. IX A[ thus altering the spec- 
trum of crystalline excitations. 

It should be emphasized that all these experiments de- 
pend crucially on the cavity's being a multimode one, and 
on the broken spatial symmetry being at least approxi- 
mately continuous. In a single-mode cavity, in which the 
broken symmetry is of the discrete (i.e., even/odd) type 
(and, moreover, the interactions are effectively infinite- 
ranged), the "solidity" is of a different kind; in particular, 
there can be neither Goldstone modes nor topological de- 
fects in the solid. 



B. Supersolid- "Mott" transition 

As discussed in, e.g., Ref. |52j, the "normal solid" state 
that competes with a supersolid is analogous to a Mott 
insulator with regard to its transport properties. In the 
setting of self-organized atom-light crystals, for laser in- 
tensities well above threshold one expects the emergent 
lattice potential to be sufficiently deep to cause the su- 
persolid state to have undergone a transition into a non- 
superfluid state. This state can be either a Mott insulat- 
ing state (which would be a normal solid) or a Bose glass 
state (which, too, would have many of the properties of 
a normal solid, particularly a non-quantized moment of 
inertia). The latter possibility arises even in the absence 
of extrinsic disorder because vacancies and dislocations 
in the self-organized lattice might dynamically generate 
disorder. The Mott insulator state (which is incompress- 
ible) and Bose glass state (which is compressible) should 
be distinguishable via, e.g., their large-scale spatial den- 
sity profiles. 

The BEC-to-Mott transition has recently been 
addressed — for the case of a single-mode cavity — in 
Ref. [14]. The case of the concentric cavity has sev- 
eral features in common with the single-mode cavity case; 
there is, however, one key difference, which follows from 
the difference in the character of the self-organization 
transitions in these two cases. Consider the phase dia- 
gram in terms of the two physically adjustable parame- 
ters, viz., the pump laser strength O (or, equivalently, the 
effective coupling constant Q and the scattering length 
of the atoms, a (which can be tuned, e.g., by approaching 
a Feshbach resonance). In the case of a single-mode cav- 
ity, the self-organization transition is continuous; there- 
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FIG. 11. Schematic zero-temperature (i.e., quantum) phase 
diagram for a BEG in a concentric cavity, with the control pa- 
rameters being the atomic scattering length a and the inverse 
effective atom-cavity coupling (^~^ (or equivalently the inverse 
laser intensity fl"^). For weak, repulsive interactions, the su- 
perfluid first undergoes seff-organization via the Brazovskii 
transition, thus forming a supersolid. If the laser intensity is 
increased further, the supersolid undergoes a transition into 
a normal solid (i.e., a Mott insulator). However, for strong, 
repulsive interactions, the uniform BEG can lose phase coher- 
ence concurrently with the first-order self-organization tran- 
sition. This situation is to be contrasted with that for the 
case of a single-mode cavity (inset), in which there should al- 
ways be a supersolid (SS) region separating the uniform fluid 
(SF) and normal solid (S) regions. First- and second-order 
transitions are marked (1) and (2) respectively. 



fore, by tuning the laser to sufficiently near threshold, 
the emergent lattice depth can be made arbitrarily small. 
Accordingly, regardless of how large the scattering length 
might be, there is always a region in which the self- 
organized lattice is too shallow to support a Mott insula- 
tor. Put differently, there is always a region of the super- 
solid phase between the liquid (i.e., the uniform BEC) 
and the normal solid. By contrast, self-organization in a 
concentric cavity occurs by means of a first-order, Bra- 
zovskii transition. The emergent lattice depth therefore 
jumps discontinuously to some nonzero value at the self- 
organization transition; if this minimum lattice depth is 
greater than that required to support a Mott insulator, 
it is possible to have a direct liquid-to-Mott transition, 
without an intermediate supersolid phase |56j . The phase 
structure of our atom-cavity system that results from the 
foregoing considerations is summarized in the schematic 
phase diagram shown in Fig. [TT] Curiously, it has the 
same morphology as that predicted for the corresponding 
liquid-to-solid transitions of *He [57] (but with temper- 
ature and pressure serving as the control parameters in 
the ^He case). 
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XI. EXPERIMENTAL FEASIBILITY 

In this section we discuss typical parameter values 
for which the Brazovskii transition, and the resulting 
self-organized state, should be observable in a mul- 
timode cavity. Recently, the self-organization transi- 
tion was observed in a single-mode cavity using Bose- 
Einstein-condensed ^^Rb atoms. The cavity was char- 
acterized by by the cavity QED parameters {g,K,j) — 
27r X (10,1,3) MHz; and the pump laser was detuned 
from the atomic resonance by A^i w 27r x 10^'^ Hz W\ . 
An appreciable photon population was observed in the 
cavity for A^ up to 2tt x 40 MHz. For these parame- 
ters, the mean-field value of the self-organization thresh- 
old r^nif is approximately 27r x 2 GHz. For a concen- 
tric cavity having these parameters, the range in fre- 
quency space of the fluctuation-dominated regime for 
the classical self-organization transition is approximately 
27r X (2 — 20) MHz. (The precise value depends on the 
healing-length parameter X: but only as x^^ i hence 
these results are not very sensitive to the choice of x-) 
As previously noted [T7], this range is greater than the 
frequency width associated either with spontaneous de- 
cay or with the intensity-noise of a typical laser (con- 
verted to frequency units). For the quantum transition, 
the fluctuation-dominated regime (for *^Rb at unit fill- 
ing far from a Feshbach resonance) would be of order 
0.1 — 1 MHz; the width of this regime can, however, be 
extended by tuning the interparticlc interaction through 
a Feshbach resonance. 

The nucleation rates, both thermal and quantal, can 
be expressed — for the case of the interacting system — as 
f/exp[— .^(fith — ^th)/i^ ~ ^th)]- In this expression, fi 
is the pump laser strength and .^ is a number of order 
unity; the expression is only valid when the exponent ex- 
ceeds unity, i.e., for 17 sufficiently near f^th- For U far 
from a Feshbach resonance, this expression implies that, 
for the regime discussed in this work, in which nucle- 
ation proceeds via the formation of large, well-defined 
droplets, the average nucleation timescale should typi- 
cally exceed the lifetime of the experiment (discussed in 
the next paragraph). It might therefore prove necessary 
to enhance V by means of a Feshbach resonance in or- 
der to explore the physics of nucleation. Note that this 
does not imply that the self-organized state is inaccessi- 
ble away from a Feshbach resonance: rather, one expects 
that the laser strength would have to be increased (de- 
creased) well past threshold before the system entered 
(left) the self-organized state. In other words, the process 
of self-organization should exhibit significant hysteresis. 

For the experimental parameters just given, Baumann 
et al. [Mj found that the lifetime of the self-organized 
BEG (i.e., the timescale on which atom loss destroys the 
BEG) was approximately 10 ms. For detunings Ac > k, 
atom loss is largely due to spontaneous scattering, which 
occurs at a rate R^ proportional to 7^7^/ A^. The corre- 
sponding timescale is long enough to enable the observa- 
tion of the family of phenomena discussed in the present 



work. In order to perform a similar experiment with cav- 
ities in the weafc-coupling regime (i.e., for smaller g), one 
would have to ensure that the spontaneous decay rate 
R~f stays below 1 kHz. This would involve satisfying two 
conditions: g'^N/Ac > 10"^ and k <^ Ac- One could do 
so by increasing the number of atoms in the BEG or by 
increasing the finesse of the cavity mirrors (which would 
reduce k). 



XII. SYSTEMS OF COUPLED LAYERS AND 
THE ORIGINS OF FRUSTRATION 

We have discussed how an equilibrium atomic cloud, 
confined by the pump laser to a plane near the equato- 
rial plane of the cavity, spontaneously crystallizes glob- 
ally into one of a family of degenerate quasi-checkerboard 
arrangements. Now let us consider an atomic cloud con- 
fined to a single plane away from the equator of the cav- 
ity. In this case, spontaneous crystallization still occurs 
but, as we shall now explain, the particular checkerboard 
arrangement into which the atoms crystallize varies stat- 
ically across the plane — energetics demands, e.g., that 
the center and edge of the cloud crystallize in distinct 
arrangements. This is a consequence of /rwstraf ion: satis- 
fying local energetic preferences introduces "fault zones" 
between locally ordered regions. 

In our analysis of equatorial-plane atomic distribu- 
tions (see Fig. Isl), we were able to focus on the family 
of degenerate modes TEM;„j having 1 = 0. To gener- 
alize our analysis beyond the equatorial plane, we must 
consider all modes that meet the degeneracy condition 
I -\- m, -\- n = KqR. Gonsider the situation illustrated 
in Fig. [I2^ , first focusing on the non-equatorial sheet 
marked (i). Near the center of this sheet, crystalliza- 
tion into modes with Z = 1 is suppressed because such 
modes have low intensity, whereas crystallization into 
I = 2 modes is favored because they have maximal in- 
tensity; away from the center, the opposite is true. The 
change in I forces a change in m or n, owing to the de- 
generacy condition, so the mode functions in the sheet 
must change across an interfacial zone between the I = 1 
and / = 2 favored regions. Therefore, either a disloca- 
tion, associated with a change in m,, or an abrupt change 
in lattice periodicity (i.e., a discommensuration) , asso- 
ciated with a change in n, is expected. (This picture 
assumes that, as is always the case near threshold, the 
self-organized lattice is not strong enough to trap the en- 
tire atomic distribution at the center or the edge of the 
sheet. The kinetic energy cost of localization, as well as 
the cost in repulsive energy, will act to spread the atomic 
cloud out.) 

Now consider a situation in which two symmetrically 
disposed sheets on opposite sides of the equator are pop- 
ulated with atoms, as in sheets (i) and (ii) of Fig. 12 
Atoms in sheet (i) and those in sheet (ii) are coupled via 
the cavity modes. Because / = 2 (^ = 1) mode func- 
tions are symmetric (antisymmetric) about the equato- 



29 




TEM 



FIG. 12. Schematic illustration of the implications of frustra- 
tion. Atoms are loaded into sheets (i) and (ii), shown as thick 
lines in panel (a), which are an integer number of pump-laser 
wavelengths apart. The dashed and dashed-dotted curves are, 
respectively, antinodal regions of the modes TEMim, which 
have low intensity near the centers of sheets (i) and (ii), and 
TEM2m, which have low intensity away from the centers of 
sheets (i) and (ii). Near the center of each sheet, atoms crys- 
tallize into the TEM2m modes; away from the center, they 
crystallize into the TEMi^' modes. Within a sheet, regions 
may be separated by faults in the ordering, as illustrated in 
panel (b). For example, on the left side the fault has the 
form of a discommensuration (see Sec. |XIl |. By contrast, the 
fault on the right side is a grain boundary. Between layers, 
the opposing parity of adjacent modes leads to frustration, 
which precludes ordering, as in the regions indicated by a A 
and a D. Grain boundaries (denoted by D) are more local- 
ized faults, and are therefore less costly, energetically, than 
discommensurations (A). 



rial plane, the atoms in the I — 2 (I — 1) arrangement 
in sheet (ii) occupy the same (opposite) checkerboard as 
those in sheet (i). If there are no dislocations, atoms 
in the interfacial zone remain disordered, because it is 
impossible for the atoms to satisfy both desiderata [or, 
equivalently, because the corresponding cavity modes in- 
terfere destructively in sheet (i) and constructively in 
sheet (ii)]. The introduction of dislocations enables the 
system to order in part of the interfacial zone, as shown 
in the right hand side of Fig. \T2}p , and is therefore pre- 
ferred. 

The full many-layer, many-mode system is expected 
to experience the same kinds of disordering effects as 
the idealization sketched above: i.e., one expects sys- 
tems slightly above threshold to develop locally crys- 
talline phases separated by zones riddled with faults. 



XIII. CONCLUDING REMARKS 

In this work, we have shown that the self-organization 
of BECs in multimode cavities is accompanied by a range 
of effects, such as fluctuation-driven nonequilibrium first- 
order transitions (both classical and quantal) , topological 
defects, rigidity, frustration, and supersolidity. We have 
developed both a nonequilibrium formalism for exploring 
such systems in general, and an effective equilibrium de- 
scription valid in the regime of greatest interest, viz., the 
quantum phase transition undergone by a Bose-Einstein 
condensate. We have outlined, moreover, how these for- 
malisms may be used to compute the correlation func- 
tions of the photons emitted from the cavity, as well as 
those of the atoms, and how such correlations may be 
detected experimentally. Finally, we have suggested re- 
alistic values of experimental parameters that might be 
used to realize self-organization and its attendant phe- 
nomenology in the laboratory. 

A major outstanding question that we hope to ad- 
dress in future work is whether the self-organized state 
in the layered three-dimensional geometry described in 
Sec. |XII| exhibits glassy dynamics. Furthermore, we see 
three promising avenues for future work extending the 
ideas developed in this paper. The first is to consider 
pumping with light sources that exhibit either thermal 
or quantum noise, instead of considering classical laser 
light. Such light sources are likely to lead to nontrivial 
effects because of the interaction between the input noise 
and the environmental noise due to photon loss. The sec- 
ond avenue is to apply the current formalism to related 
instabilities of atoms in cavities: e.g., the phenomena 
of collective atomic recoil lasing (CARL) |58j and "ex- 
cess noise" [55]. Such phenomena are well understood 
at the mean-field level, but the present formalism makes 
it straightforward to address the effects of fluctuations. 
Finally, it might be of interest to explore the physics of 
many fermions trapped in multimode cavities, especially 
the prospects for their Cooper pairing. There are two 
possible ways to arrange this: the first is to exploit the 
previously mentioned analogy (see Sec. Ill I between pho- 
tons and phonons in order to achieve photon-mediated 
Cooper pairing; the second is to use the phonons of the 
emergent crystalline phase as the "glue" that would bind 
the Cooper pairs. 
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In terms of Ai and A2, the Hamiltonian H has the form 



Appendix A: Toy equilibrium model 

In this Appendix we elaborate on the connection be- 
tween quantum phase transitions in quasi-equihbrium 
cavity QED settings and the conventional theory of quan- 
tum phase transitions as level-crossings (as discussed, 
e.g., in Ref. f5(T), by means of a toy model. This model 
has the following components: (i) a laser mode a^ of 
frequency oj^, (ii) a cavity mode ac of wavevector k and 
frequency toc = i-ol + ^c, and (iii) an anharmonic atomic 
phonon b of wavevector k and frequency cOp . The Hamil- 
tonian for the model is given by: 



- 2 



A^ + 472x2 

(A5) 

Expanding in powers of x, one finds that the coefficient 
of x^ in the Hamiltonian is now given by 



1 + (4^1-4^2)1^. 

Ac 



(A6) 



For a fixed N , the smallest value that the last term can 
attain is 



H = hujj^aj^a^ 






hutpb'b 



hT{ala^+h.c.){b'^ 



^Mb^bf 

b), (Al) 



where F is the photon-phonon coupling constant, and 
A is a parameter describing the strength of the an- 
harmonicity of the phonon. In what follows we shall 
rescale all energies by hwp in order to de-dimensionalize 
them; the dimensionless parameters will be denoted as 
lJl etc. As H commutes with the total number of pho- 
tons A/'(= aj^fli + al^ac), one can diagonalize it in a space 
of fixed total photon number Af. (For sufficiently large 
A^, the difference between number and coherent states is 
irrelevant.) Let us attempt to find the ground state for 
a particular value of Af. It is helpful to simplify further 
and treat the phonons as being "classical" by taking the 
commutator [6^,6] = 0; this is equivalent to rewriting 
b — ^/mujj2h[x + i{p/rnu!)] and taking the to — >■ 00 limit. 
In this limit, all terms depending on p in the Hamiltonian 
are suppressed; thus the Hamiltonian can be rewritten, 
in terms of a; oc {b + b^), as 



H = Wi(4aL + ac«L) + ^L ' 



+x^ + Xx^ + 2'yx{al;aj^+h.c.). (A2) 



,t 



In this expression the prefactors accompanying x have 
been absorbed as appropriate into the various coupling 
constants. The first term is simply a constant times Af, 
and can be ignored for our purposes. As x is also now a 
good quantum number, it is now possible to diagonalize 
H in a manifold of fixed Af and x. This can be achieved 
via canonical transformation from from aj^/c to Ai/2'. 



fli == aAi + v 1 ^ ct^A2, 
ac = — V 1 — 0:2^1 + aA2, 



where 



1± 



^C 



472^2 + A^ 



(A3a) 
(A3b) 



(A4) 



Ac 



(A7) 



Thus, when Gj'^/Ac < l/A/", the coefficient of x^ is al- 
ways positive; hence, the ground state always lies in the 
X — manifold. However, for 6j'^/Ac > l/A/", the 
ground state is that in which AIA2 = Af, and 



Ac 



2A 



(A8) 



Hence, as 7 is increased, the ground-state value of x (for 
fixed Af) becomes nonzero, singularly. In other words, 
the lowest level of the x = manifold and that of an 
X 7^ manifold cross at the critical value of 7: such 
a level crossing is known as a "quantum phase transi- 
tion" [5D] . Although the phenomenon of self-organization 
in a multimode cavity involves substantially more than a 
single photon and a single phonon, its essential character 
in the equilibrium limit, for a classical laser, is also that 
of a ground-state level-crossing occurring at some fixed 
photon number Af. 



Appendix B: Determining the healing length 
parameter x 

The parameter 1/x, which sets the healing length for 
crystallinity, is a measure of how weak the atomic cou- 
pling (imn to I ^ modes is, relative to the coupling 
to Z = modes. The value of x is determined by the 
following effects: (1) the atoms, being confined near the 
equatorial plane of the cavity, couple most strongly to 
modes that have the highest amplitude there, and thus 
to the lowest-order modes along the z direction (i.e., g is 
in effect a function of I once S(x) is projected onto the 
equatorial plane); (2) higher-order modes along the z di- 
rection have lower finesse and hence couple more weakly 
to the atoms (i.e., k is a function of I); (3) the effec- 
tive laser-cavity detuning, Ac — g^N/ Aa, is a function 
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of / because g is; and (4) for a nearly concentric (or con- 
focal, planar, etc.) cavity — the experimentally relevant 
case — Ac is a function of I because higher-order modes 
have lower resonant frequencies. Thus, C is in general a 
complicated function of /. 

For the specific case of a concentric cavity, one can gen- 
eralize the mode structure given in Sec. |VIl"5] by impos- 
ing, e.g., finite- well rather than hard- wall boundary con- 
ditions on 9. In this case, it is clear that both 5;/.go ^nd 
hi/kq (in which the subscripts denote the appropriate 
value oil) must go approximately as (1— constx/^), where 
the constant is approximately 1/^maxi i-^-; the number of 
higher-order modes that are of sufficiently high finesse 
to couple significantly to the atoms. Thus, x can be 
taken to be approximately 1/^maxi ^P to a factor of or- 
der unity; the physically relevant quantities ^o and flth, 
which depend on x^/"^ and x^^^ respectively, should not 
be sensitive to this neglected factor. 

Note that this analysis ignores effect (4). According 
to Eq. (19.24) of Ref. [3D], this term would lead to a fre- 
quency shift Ac(/)/Ac(0) ex —I, with a proportionality 
constant depending on the distance from concentricity. 
Thus, C should have a positive linear term in /, in addi- 
tion to the negative quadratic term; the potential effect 
of such a term is to favor some family oi I = Iq y^ modes, 
even at the equator, over I = ones. Similar considera- 
tions apply to other geometries, such as the confocal or 
planar cavities. 



Appendix C: Renormalization-group flow equations 
for the quantum Brazovskii model 

In this Appendix we outline our derivation of the 
renormalization-group (RG) flow equations for the quan- 
tum Brazovskii model discussed in Sec. |VIIF Our pro- 
cedure parallels that discussed in App. A of Ref. [29] for 
the classical Brazovskii model; that work, in turn, was 
based on techniques developed by Shankar ^21j for the 
Fermi liquid. The objective of the renormalization-group 
procedure is to arrive at a spatially coarse-grained effec- 
tive theory in terms of modes for which m + n ~ Aq. 
This is done by progressively integrating out modes for 



which \m 



An I > B, where B is referred to 



the renormalization-group "scale" The microscopic, or 
"bare" theory has an RG scale that is associated with 
the physical high-energy cutoff (e.g., A^ in the cavity 
QED case); this scale is progressively decreased by the 
integrating out of "shells" of modes, i.e., modes for which 
-Sncw < |m -I- n — AqI < -Bold, thus yielding effective the- 
ories involving progressively fewer modes. As one inte- 
grates out these shells of modes, the theory maintains its 
basic structure, but the various coupling constants flow; 
therefore, the coupling constants are functions of B. For 
example, the microscopic values of the parameters, TZ 
and W, are their values at a value of B determined by 
A^, whereas the fully coarse-grained parameters r and 
u are those corresponding to i? = 0. At each step in the 





FIG. 13. Feynman diagrams involving the six-point vertices 
associated with the couplings Wi {i — 1,2,3), to one-loop 
order. The six-point vertices are denoted as grey squares. 
The diagram shown on the left renormalizes the four-point 
vertices; that on the right renormalizes the six-point vertices. 



RG procedure, we integrate out all the frequency compo- 
nents associated with that spatial shell; in making this 
choice [which has the advantage of preserving the causal- 
ity structure of the microscopic theory (as mentioned in 
Ref. [S])] we are following Refs. [HJBl]. 

We begin with the observation that the action can be 
written as follows: 



S=S2+S4 + ' 



(CI) 



where 



52 = ^y'dlr(l)$(l)$(l), (G2) 

Si = ^^ /"dld2rf3d4M(1234)$(l)$(2)$(3)$(4) 
x5(l + 2 + 3 + 4), 

and notation of the form dl is to be interpreted in the 
following way: 



dl = 



(2^) 



A^drj 



|m+n-Ao|<S 



duj / de 



(C3) 

where 77 and 9 are, respectively, the quasi-radial and 
quasi-angular variables discussed in Sec. |VIIG and are 
treated here as continuous. We now perform the RG 
transformation, which involves two steps: (i) integrating 
out all modes satisfying B/b < |to -I- n — Ao| < B for 
6=1-1-/ with / <C 1; and (ii) rescaling the spatial coordi- 
nates and the flelds in the action so that the new action 
is similar in form to its predecessors under the RG trans- 
formation. Under this rescaling, 77 -> 677 and $ -^ ^/b. 
Under the combined effects of integrating out the shell 
of modes (at one-loop order, via the diagram in Fig. [3]) 
and the rescaling, the coefficient of the quadratic term 
transforms as follows: 



r{B) ^ r{B/b) = b'(r{B) + A2) 



(G4) 



where 
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A2 = aui{B) / d-q—-— 



(C5) 



By setting b = l + l and differentiating with respect to /, 
we arrive at a differential version of tlie RG equation for 
r, viz., 






B2 



(C6) 



We have introduced the notation ui instead of u because, 



as we shall now see (and as we anticipated in Sec. VII D I, 



the renornialized value of u depends on whether the four 
9s entering the vertex are identical or not. The generic 
vertex, in which the four values of 9 are not all identical, 
is associated with the coupling ui] the special vertex, in 
which the four values of 9 are all identical, is associated 
with the coupling U2. Similarly, we shall denote the six- 
point vertices respectively as wi, W2, and w^, depending 
on whether two, four, or all six of the incoming propa- 
gators have identical values of 9. The RG equations for 
(ui,'U2) can be derived via the same procedure as those 
for r; they read as follows: 



is as follows: 

dr 

dB 

dui 

~dB 

dU2 

~dB 
dwi 
~dB 

dw2 
~dB 

dvh 
dB 



Ml 



(f + e2)l/2 



Wl 



^ + 62 

2ul 



(f 



-62)1/2' 
W2 



^2)1/2' 

SuiWi 

(f 4-i52)3/2 + f + ^2 ' 
U2W2 



^2 (? 

2u? 



4u? 



(f-hS2)3/2 

12uf 



+ 



8uiW2 



AuiWi 
r + B^' 



(f + 62)3/2 f + B^' 



(C13) 
(CM) 
(C15) 
(C16) 

(C17) 
(C18) 



The appropriate microscopic values (i.e., initial 
conditions for these equations) are as follows: 
(f(cx)),ui 2(00), Wi_2, 3(00)) = (7?,, 1,0). Numerically 
integrating the equations yields the coarse-grained 
parameters, which are plotted in Fig. [5] and discussed in 
the main text. 



Appendix D: Effective temperatures 

In this Appendix we briefly explain why, and for what 
purposes, k in Sec. |VIIl| behaves as an effective tempera- 
ture. The basic result we shall review is as follows: con- 
sider a quadratic Keldysh action that has a q — q compo- 
nent of the (low-frequency) form 



dui 

~dr 

du2 



2ui 



= 2^9-2 



aufB 
r{B) + B^ 

au\B 
r{B) + B^ 



awiB 



Vr{B)+B^' 
aw2B 
^r{B)+B^' 



(C7) 
(C8) 



Similar equations can be derived for the three w cou- 
plings; these involve the diagrams in Fig.|3ji and Fig. 13 
Before writing these down, however, we introduce the 
following convenient changes of variables, which serve to 
de-dimensionalize the RG equations (following Ref. [2 9 j ) : 



i^r^, / duj (j)^^q{uj) (t)^,^, 



-UJ 



(Dl) 



for some real set of parameters F. (The generic set of 
indices parameterized by v can describe positions, mo- 
menta, mode indices, etc.) The long-time dynamics of 
such a theory can then be described by a Langevin equa- 
tion, having a white noise term of strength F. 

In what follows, we shall suppress the ly index; the 
argument, which is adapted from Ref. |31) . can be made 
independently for each v. First, note that Eq. (Dl) is 
expressed in the time domain as 



B = B{aU) 



-1/2, 



r{B) = r{l){aU)-^ e-^\ 
Mi,2(S) = Mi,2(0We-2', 

Wi^2,z{B) = 1^1,2,3(0^"^ "^e 



(C9) 
(CIO) 
(Cll) 
(C12) 



iV / duj(j)q{t)(l)q{t). 



(D2) 



This term appears in the system's partition function Z 
J D(j)q{t) D(f>c{t) expiS, in the form 



Z = / Dq^q{t) exp -F / dtcbqit) 4>q{t) X 



(D3) 



in which the ellipses denote factors resulting from other 
terms in the action. Eq. (D3) above can be rewritten, by 



In terms of the new variables, the full set of RG equations means of a Hubbard-Stratonovich transformation, as 
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j D<j,q{t) D^{t) e^^ 



m' 



(D4) 
Once this is done, the full partition function, which also 
includes terms linear in (pq (from the retarded and ad- 
vanced components of the Keldysh action) can be written 
as follows: 



i?C(t)D(/.,(i)exp 



dtm' (D6) 



x5[Ut)G{t,t')~at')]- 

Accordingly, the dynamics of the system is described, at 
long times, by a sum over classical histories in the pres- 
ence of a Langevin white-noise term ^, the fluctuations 
of which are given by 



{mat'))^v5{t-t'). 



(D7) 



Z = I Dat)DMt)c^p(-f J dtatf] (D5) 



where G represents some (unspecified) integral kernel 
that couples (/>c and (j^q- Note that the c — c term in the 
action is absent, via causality, as discussed in Ref. [ST] , 
If one now integrates out 0g, one finds that the partition 
function is given by 



It follows that the long-time dynamics of the system is 
classical rather than quantal; thus any phase transition 
that the system undergoes is a thermal rather than a 
quantum phase transition. By contrast, for any system 
that undergoes a true quantum phase transition, the co- 
efhcient of the q — q component of the Keldysh action 
vanishes at low frequencies, typically as a power law, 
IcliJ" (see, e.g., Ref. !44J. This corresponds to power-law 
decay of noise correlations in the time domain. 

A closely analogous argument (see, e.g., Ref. [3T] ) 
shows that, assuming the q — q component of the action is 
frequency-independent as a; — >■ 0, the rate of escape from 
a metastable state is given by an effective Arrhenius for- 
mula with a temperature proportional to F. 
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